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ABSTRACT 

Using deep near-infrared and optical observations of the Hubble Deep Field North from the Hubble 
Space Telescope NICMOS and WFPC2 instruments and from the ground, we examine the spectral energy 
distributions of Lyman break galaxies (LBGs) at 2.0 ^ z < 3.5 in order to investigate their stellar 
population properties. The ultraviolet-to-optical rest-frame spectral energy distributions (SEDs) of the 
galaxies are much bluer than those of present-day spiral and elliptical galaxies, and are generally similar 
to those of local starburst galaxies with modest amounts of reddening. We use stellar population synthesis 
models to study the properties of the stars that dominate the light from LBGs. Under the assumption 
that the star-formation rate is continuous or decreasing with time, the best-fitting models provide a lower 
bound on the LBG mass estimates. LBGs with "i*" UV luminosities are estimated to have minimum 
stellar masses ~ lO^*' -Mq, or roughly 1/lOth that of a present-day L* galaxy, similar to the mass of 
the Milky Way bulge. By considering the photometric effects of a second stellar-population component 
of maximally-old stars, we set an upper bound on the stellar masses that is ~ 3 — 8x the minimum 
mass estimate. The stellar masses derived for bright LBGs are similar to published estimates of their 
dynamical masses based on nebular-emission-line widths, suggesting that such kinematic measurements 
may substantially underestimate the total masses of the dark matter halos. We find only loose constraints 
on the individual galaxy ages, extinction, metallicities, initial mass functions, and prior star-formation 
histories. Most LBGs are well fit by models with population ages that range from 30 Myr to ~ 1 Gyr, 
although for models with sub-solar metallicities a significant minority of galaxies are well fit by very 
young (< 10 Myr), very dusty stellar populations, A(1700 A)> 2.5 mag. We find no galaxies whose SEDs 
are consistent with young (< 10* yr), dust-free objects, which suggests that LBGs are not dominated by 
"first generation" stars, and that such objects are rare at these redshifts. We also find that the typical 
ages for the observed star-formation events are significantly younger than the time interval covered by 
this redshift range 1.5 Gyr). From this, and from the relative absence of candidates for quiescent, 
non-star-forming galaxies at these redshifts in the NICMOS data that might correspond to the fading 
remnants of galaxies formed at higher redshift, we suggest that star formation in LBGs may be recurrent, 
with short duty cycles and a timescale between star-formation events of < 1 Gyr. 

Subject headings: early universe — galaxies: evolution — galaxies: stellar content — galaxies: 
starburst — infrared: galaxies 



1. INTRODUCTION 

The last few years have seen rapid advances in the study 
of galaxies at very high redshifts, especially aX z > 2. In a 
large part, this has been due to the development of simple 
photometric techniques to select high redshift objects by 
characteristic color signatures. Neutral hydrogen absorp- 
tion within galaxies and from the intergalactic medium 
strongly attenuates flux shortward of Lyman a (1216 A) 
and the 912 A Lyman limit, producing a spectral "break" 
that provides a color signature to identify objects at high 
redshift. Such color selection methods have long been used 
to identify high redshift QSOs (e.g., Warren et al. 1987), 
and were latter applied to set limits on the number of star- 



forming, faint galaxies at z « 3 (Guhathakurta, Tyson, & 
Majewski 1990; Songaha, Cowie, & Lilly 1990). Using this 
method with a custom suite of broad-band filters, Stei- 
del & Hamilton (1992) and Steidel, Pettini, & Hamilton 
(1995) established the existence of a significant number of 
high redshift galaxy candidates, now commonly described 
as Lyman break galaxies (LBGs), which were then spectro- 
scopically confirmed to have z « 3 (Steidel et al. 1996a, b). 
By now, nearly 1000 LBGs have measured redshifts, pro- 
viding a rich sample for studying the properties of star- 
forming galaxies at high redshift. 

Understanding the epoch at which galaxies assembled 
the bulk of their stellar mass remains one of the key cos- 
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mological questions. Although new surveys are produe- 
ing a growing inventory of photometric and spectroseopie 
data on LBGs, it is not yet clear how these objects fit 
into the ancestral history of the present-day galaxy pop- 
ulation. The process of galaxy assembly depends on the 
particular parameters of the cosmological model, e.g., the 
biasing of baryonic mass relative to the underlying dark 
matter distribution, and the physics and feedback of star 
formation. In hierarchical models, galaxy formation is a 
continuous process driven by mergers of lower mass sub- 
components over a wide range of redshifts. In other sce- 
narios, star formation and galaxy assembly takes place on 
rapid timescales at high redshift, with the galaxies evolv- 
ing in situ thereafter, more or less passively, to the present 
day. These different models predict substantially different 
histories and timescales for the formation and assembly 
of early stellar populations, especially at ^; > 1. Ob- 
servational constraints of the mass assembly history are 
inconclusive. The presence of old stars (> 10 Gyr) in 
the Galactic bulge and the spheroidal components of M31 
(e.g., Renzini 1999, Rich & McWilliam 2000), the exis- 
tence of a population of old ellipticals at z > 1 (Dickinson 
1995; Dimlop et al. 1996; Daddi, Cimatti, & Renzini 2000), 
and the small scatter in the color-magnitude relation of 
cluster ellipticals to 2; ~ 1 (Stanford, Eisenhardt, & Dick- 
inson 1998) all support the picture of rapid collapse and 
enrichment at high redshifts, and imply that these systems 
formed a significant fraction of their stellar mass by z > 2. 
Conversely, the observed clustering properties of galaxies, 
the steep faint-end slope of the local luminosity function 
(e.g., Folkes et al. 1999, and references therein), and evi- 
dence for the relative paucity of massive galaxies at z > 1 
(Cowic et al. 199G; Kauffman & Chariot 1998) have been 
cited to support the hierarchical picture. 

One important and unresolved issue is the stellar pop- 
ulation content of LBGs, and their stellar and dark mat- 
ter masses. Stcidel and collaborators (e.g., Steidel et al. 
1996b; Pettini et al. 2000) have found that the spectra of 
LBGs are broadly consistent with ongoing star formation, 
and are similar to those of local starburst galaxies. Gi- 
avalisco, Steidel, & Macchetto (1996) and Lowenthal et al. 
(1997) have emphasized the compact sizes of LBGs, with 
half-light radii of only a few kpc. The strong clustering 
observed for LBGs favors relatively large dark matter halo 
masses (Giavalisco et al. 1998; Adelberger et al. 1998). 
These observations are consistent with simulations (e.g., 
Governato et al. 1998; Baugh et al. 1998) that predict 
that LBGs evolve to become the bulges and spheroidal 
components of present-day, high-mass (~ L*) galaxies, 
preferentially in dense environments. In contrast, Lowen- 
thal et al. (1997) argued that if z ^ 3 LBGs were to fade 
without continued star formation or subsequent mergers, 
then their sizes and luminosities would ultimately resemble 
those of present-day, low-mass dwarf elliptical/spheroidal 
galaxies. They suggest that the LBGs are possibly low- 
mass objects, imdergoing intense, brief periods of intense 
star formation, which place them above the magnitude 
limit and within the color-selection criteria. Observations 
of nebular emission-line widths for LBGs find typical val- 
ues in the range cr ~ 60 — 120 km s^^ (Pettini et al. 1998, 
2001; Teplitz et al. 2000a; Moorwood et al. 2000). Com- 
bined with the sizes for the galaxies, these line widths 



suggest virial masses ^ although it is unclear 

whether such measurements really trace the total mass of 
the dark matter halos of these galaxies, which may extend 
well beyond the limits of the star forming, UV-bright, line 
emitting region. 

Optical photometry has shown that the rest-frame ul- 
traviolet (UV) continua from LBGs are typically redder 
than is expected from ongoing, unreddened star forma- 
tion, suggesting the presence of dust (e.g., Meurer et al. 
1997; Dickinson 1998; Sawicki & Yee 1998; Meurer, Heck- 
man, & Calzetti 1999; Steidel et al. 1999; Adelberger & 
Steidel 2000). Nebular-line measurements for LBGs, men- 
tioned above, have sometimes yielded estimates of star- 
formation rates (SFRs) that are higher than those derived 
from the UV continuum measurements (but see Pettini et 
al. 2001), again suggesting some amount of dust extinction 
(A1700 = 1-2 mag). 

Most photometric studies of LBGs have focused on their 
rest-frame UV light, which is dominated by short lived, 
massive stars. Limited to this portion of the spectral en- 
ergy distribution (SED), it is not possible to disentan- 
gle effects of metallicity, age, extinction, and past star- 
formation history. By extending the photometric base- 
line to rest-frame optical wavelengths using near infrared 
(NIR) data, we may measure light from both short and 
longer-lived stars, and may hope to constrain the stellar 
population mix in these objects. For all but the youngest 
galaxies, the optical rest frame spectrum has a strong con- 
tribution from the light of later-type, main-sequence stars 
(A-type and later) and evolved red giants, and we may 
therefore expect that NIR measurements will be more sen- 
sitive to the total stellar mass, which is dominated by 
lower mass stars. Also, the longer wavelengths arc less 
affected by dust extinction that may strongly attenuate 
ultraviolet light. 

Sawicki & Yee (1998) have studied the SEDs of 17 LBGs 
from the Hubble Deep Field North (HDF-N, Wilhams et 
al. 1996) using optical WFPC2 photometry augmented 
by ground-based infrared JHKg measurements from data 
obtained by Dickinson and collaborators (cf. Dickinson 
1998). They concluded that these objects are generally 
quite young (median age ^ 25 Myr) and highly reddened 
[median extinction E{B — U) ~ 0.3, ^1700 — 3.2 mag]. 
Based on the derived SFRs and timescales, they concluded 
that the light from LBGs is dominated by short bursts of 
star formation and will likely only produce ^ 5% of the 
total stellar mass of a present-day L* galaxy. However, 
the HDF-N LBGs are typically quite faint, and the avail- 
able ground-based infrared data provide measurements 
with rather low signal-to-noise ratios (S/N). Addition- 
ally, there are special challenges to matching the ground 
based and WFPC2 photometry given the different resolu- 
tions and S/N of the images. 

In this paper, we investigate the properties of a sam- 
ple of 33 LBGs with spectroscopically confirmed redshifts 
2.0 < z < 3.5 drawn from the HDF-N, using new data 
from the Near-Infrared Camera and Multiobject Spectro- 
graph (NICMOS, Thompson et al. 1998) on board the 
Huhhle Space Telescope (HST). In §2. we present the 
broad-band photometry used for the ensemble of objects. 
In §3, we compare the LBG SEDs to empirical spectral 
templates for various galaxy types. In §4, we constrain 
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the stellar population mix in these objects by fitting the 
observed SEDs with stellar population synthesis models. 
In §5, we discuss the implications for the galactic stel- 
lar masses and galaxy formation/evolution scenarios from 
the best-fit models. Finally, in §6, we present the con- 
clusions. Throughout this paper, we use a cosmology 
with nM = 0.3, A = 0.7, and h = (iJo/100 km s'^ 
Mpc-i) = 0.7. At the typical redshift z w 2.7 for the 
galaxies studied here, masses and SFRs derived from our 
model fitting procedure would be ^ 50% smaller in an 
Einstein-de Sitter universe with the same Hubble con- 
stant, and ~ 10% smaller for an Qm = 0.3 universe with 
A = 0. We present all magnitudes in the AB system, 
rriAB = 31.4 — 2.51og(/^/l nJy). We will denote galaxy 
magnitudes from the WFPC2 and NICMOS bandpasses, 
F300W, F450W, F606W, F814W, FUOW and F160W as 
U300, -B450, ^606, -^814, Jiw, and i?i6o, respectively. 

2. THE DATA AND THE LYMAN BREAK GALAXY SAMPLE 

The original HDF-N WFPC2 observations provide 
high-quality photometry of faint galaxies in four band- 
passes from 0.3-0.8 /xm. At ^; > 2, the Lyman limit is 
shifted into the WFPC2 F300W filter, permitting color 
selection of UV bright star forming galaxies with 2 < z < 
3.5. The WFPC2 images sample rest-frame UV wave- 
lengths at these redshifts (Aq ^ 1000 — 2000 A), and are 
thus primarily sensitive to light from massive young stars 
modulated by dust extinction. To study the photometric 
properties of the LBGs at rest-frame optical wavelengths, 
where longer-lived stars may contribute to the luminosi- 
ties, where dust extinction is less severe, and where the 
most is known about the low- redshift galaxy population, 
we must observe in the near infrared. This was one of 
the primary justifications for our NICMOS survey of the 
HDF-N, in which we mapped the complete WFPC2 field 
of view at 1.1 /xm and 1.6 /zm. 

The NICMOS data reach rest frame optical wavelengths 
A > 0.4 /Lim only for galaxies with z < 3. Observations at 
still longer wavelengths are desirable for many reasons. For 
galaxies at 2 < z < 3.5, the Kg band (2.16 /im) samples 
the light from Aq ~ 5000 — 7500 A, i.e., well into the optical 
rest frame, providing a better "lever arm" to measure the 
properties of lower mass, longer-lived stars. Deep obser- 
vations with HST longward of ~ 2 /zm are impractical due 
to the warm telescope assembly. Therefore, we extended 
our data to the ifg-band using observations with the In- 
frared Imager (IRIM) at the KPNO 4m Mayall telescope 
(Dickinson 1998; Dickinson et al. 2001). Using a technique 
developed by Fernandez-Soto, Lanzetta, & Yahil (1999), 
we have analyzed the Kg image to optimally extract pho- 
tometry matched to the WFPC2 and NICMOS data. The 
details of this method will be presented elsewhere (Pa- 
povich & Dickinson 2001), but in summary, we use the 
HiQo data to create two-dimensional image templates for 
each object, convolve these to match the point spread func- 
tion (PSF) of the ground-based data, and finally scale and 
fit the convolved templates to the Ks-ha.nd image to ex- 
tract fluxes. This method eliminates concerns about PSF 
and aperture matching effects on the relative photometry 
between the HST and ground-based images, and permits 
deblending of images for objects partially merged by the 
ground-based seeing. We have performed detailed simula- 



tions to test the reliability of this technique and to under- 
stand the distribution of fiux uncertainties. In both cases 
where the objects are well detected or are below the flux 
limit of the ground-based image, the method measures ro- 
bust fluxes or upper limits. However, most of the LBGs 
used in this work are well detected at Kg. Although the 
ifg-band data are not as deep as the HST data, they pro- 
vide our only access to the rest-frame optical wavelengths 
for z >3 galaxies. 

A detailed discussion of the infrared observations and 
photometric catalogs will be presented in detail elsewhere 
(Dickinson et al. 2001, see also Dickinson 1999, 2000); a 
summary of the data used here is given in Table 1. Briefly, 
we resampled the NICMOS, WFPC2, and IRIM images 
to the same plate scale (O'.'OS pixcl^^) and convolved the 
WFPC2 and Jno images to match the PSF of the Hiqq im- 
age. Our tests indicate that this PSF convolution matches 
point source photometry between bands to within 5% for 
aperture radii > O'.'l. Object detection and photometry 
were done with SExtractor (Bertin & Arnouts 1996) on a 
combined FllOW -|- F160W image. Relative photometry 
from the HST WFPC2 and NICMOS data was measured 
through matched isophotal apertures defined on the NIC- 
MOS images. The isophotal apertures were used for the 
relative colors because our tests indicate that these most 
closely match the Kg photometry derived using by the 
template fitting method (see Papovich & Dickinson 2001). 
In order to correct for galaxy light outside the isopho- 
tal apertures, we then scaled the isophotal photometry 
for each object by the ratio of the flux within an ellipti- 
cal aperture defined by radial moments of each galaxies' 
F160W light profile (i.e., the SExtractor "MAG_AUTO'' 
measurements) to the isophotal flux. The flux correction 
was < 20% for nearly all of the galaxies, with a median 
value of 4.4%. 

In order to ensure that the NICMOS and Kg photom- 
etry reach optical rest frame wavelengths, we limit our 
galaxy sample to a redshift interval similar to that selected 
by the "C/300-dropout" Lyman break criteria. Extensive 
spectroscopy from the Keck telescope has measured red- 
shifts for a subset of HDF N LBGs (mostly from Stcidel 
et al. 1996, Lowenthal et al. 1997, Dickinson 1998, and 
Cohen et al. 2000). In order to remove the redshift as a 
variable when fitting and analyzing the spectral energy dis- 
tributions of the galaxies, we restrict our primary analysis 
to those galaxies with spectroscopic redshifts in the range 
2-0 ^ z < 3.5. These objects are listed in Table 2. In 
several cases, our catalogs have split objects into separate 
entries that might be regarded as a single object with com- 
plex, multiple structure. For one object (HDF 2-239, us- 
ing the catalog number from Williams et al. 1996) we have 
merged two very faint components with the brighter, main 
body of the galaxy. In general, however, we have left the 
pieces separate. In several cases, the pieces have notably 
different colors, and thus independent analysis may be in- 
structive. In §4.3 we note that two galaxies have SEDs 
that seem inconsistent with their reported spectroscopic 
redshifts, and thus disregard them from the fitting analy- 
sis. We also exclude one other object, HDF 4-852.12 (NIC 
824), which is apparently a broad-line AGN at z — 3.479 
(Cohen 2001). 

Because all the objects in our subsample are bright 
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Fig. 1. — Color-magnitude diagram for galaxies at 2 < 2 < 3.5 in the HDF-N. The redshift distribution has been divided into two 
subsamples: 1.9 < z < 2.7 (blue circles), and 2.7 < z < 3.5 (red triangles). Filled symbols show the spectroscopically confirmed LBGs in 
our sample, while open symbols are HDF-N galaxies with photometric redshifts 1.95 < Zphot < 3.50 (Budavari et al. 2000; Dickinson et al. 
2001). The line shows the color— magnitude evolution as a function of age (labeled) for a JVl = 10^*^ -Mq stellar population with a 5-function 
star— formation history, and observed at 2 = 2.7. There are few photometric candidates for red, non-star-forming HDF galaxies at these 
redshifts. One possible candidate, the "J— dropout" object HDF-N J123656. 3+621322 is indicated. The S/N for the Vgoe measurement on 
this object is < 2, and thus the Vgoe ~ ^^^160 color is better considered as a lower limit. Although its nature and redshift are unknown, this 
object might plausibly be an example of a non-star— forming or dust— obscured galaxy in this redshift range (Dickinson et al. 2000). 



enough for spectroscopy, they represent only the bright 
end of the luminosity function. They therefore may be 
a biased subsample, either by color (i.e., the most UV- 
bright and actively star-forming galaxies) or possibly by 
mass. However, it is important to remember that the 
HDF-N LBGs are on average less luminous than typ- 
ical LBGs in the ground-based samples of Steidel and 
collaborators. In Figure 1, we show a color-magnitude 
diagram for all HDF-N galaxies selected from the NIC- 
MOS F110W+F160W data with Hiqq < 26.5 and with 
either spectroscopic or photometric redshifts (the latter 
from Budavari et al. 2000, using the NICMOS and Ks 
infrared photometry) in the range 1.95 < z < 3.5. In 
principle, the infrared-selected photometric redshift sam- 
ple could identify red galaxies with little or no UV flux 
that might be missed in optically selected Lyman break 
samples. However, as noted by Dickinson (2000), there 
are actually very few photometric candidates for galaxies 
in this redshift range that would not be identified in an op- 
tically selected sample. The color-magnitude distribution 
shows most HDF-N galaxies in this redshift range have 
broadly similar Vgoe ~ ^^160 colors, and that the spectro- 
scopically confirmed LBGs span the range of Vgoe — Hiqq 
colors observed for the large majority of objects in the 
photometric redshift sample. Figure 1 includes a line in- 
dicating the colors and magnitudes expected as a function 
of age ai z — 2.7 for a stellar population formed in a single 
burst with a stellar mass 10^" Mq- There arc a few galax- 
ies Veoe ^ Hieo colors slightly redder than the majority of 



the LBGs (but which are still well detected in the WFPC2 
data), including interesting objects such as the radio and 
X-ray source J123651. 7+621221 (Richards et al. 1998, 
Dickinson et al. 2000, Hornschemeier et al. 2000, Brandt et 
al. 2001), with photometric redshift Zphot ~ 2.6. One very 
red object, the so-called "J-dropout" J123656. 3-1-621322, 
is only marginally detected in the WFPC2 images and 
might plausibly be a reddened or non-star-forming galaxy 
at z > 2, or more speculatively an object at z > 10 (Dick- 
inson et al. 2000). However, with this possible exception, 
of the « 200 galaxies with (photometric or spectroscopic) 
redshifts, 2.0 < z < 3.5, there are no other good candi- 
dates for relatively massive, non-star-forming HDF galax- 
ies at z > 2 and with ages > 0.3 Gyr. 

We note a color-magnitude trend in Figure 1, in the 
sense that the brightest galaxies at Hiqq tend to have red- 
der optical-infrared colors. Because this is an infrared- 
selected, complete sample, there is no reason to think 
that this is due to any selection effect. To test the ro- 
bustness of this correlation, we performed a Monte Carlo 
boot-strap test by randomly reassigning colors from the 
^606 ~ HiQQ distribution (with substitution) to the -ffieo 
magnitudes. We find that the observed trend is randomly 
reproduced with a linear slope as steep or steeper in only 
~ 0.03% of the Monte Carlo simulations, which supports 
a strong likelihood that this correlation is not obtained 
from random sampling. This trend might suggest that the 
most luminous galaxies in this redshift range are compara- 
tively "older" (i.e., have lower ongoing SFRs for their total 
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Fig. 2. — Ensemble photometry of the 33 LBGs compared to empirical galaxy templates. The green data points are the measured LBG fluxes 
from the WFPC2 iU3oo,Bi5o,V606,hl4), NICMOS { Jiio,-ffl6o), and IRIM (Ks) images, normalized to common flux density at F^(2000 A) 
(black triangles indicate upper limits). Top panel: The yellow shaded region is defined by spectral templates for ordinary galaxies from 
Coleman et al. (1980), spanning the range from the earliest types (elliptical), through Sab, Sbc, and Scd spirals, to the bluest, Magellanic 
Irregular (Im) templates. While the Im template is similar to the average LBG, the earlier— type, Hubble sequence galaxies are all much 
redder than any of the LBGs. Bottom Panel: The yellow shaded region spans the range of empirical starburst galaxy SEDs from Kinney et al. 
(1996), from the bluest [NGC 1705, E{B - V) ~ 0.0] to reddest [0.6 < E(B - V) < 0.7] templates. Common spectral features arc indicated 
by dotted red lines. The LBG SEDs fit comfortably within the envelope defined by the starburst templates. 



mass), or that they are dustier. A similar color-magnitude 
trend has been noted for the UV light from HDF LBGs 
(IMeurer et al. 1999), and interpreted as a measure of vary- 
ing dust content and/or metallicity. 

3. COMPARISONS TO EMPIRICAL GALAXY SPECTRA 

The HDF-N Lyman break galaxies have significantly 
bluer colors than those of local spiral and elliptical galax- 
ies. In Figure 2, we compare the SEDs of the HDF-N 
LBGs with empirical UV-to-optical spectral templates 
for local galaxies from Coleman, Weedman, & Wu 1980 
(henceforth CWW). The LBG SEDs are uniformly bluer 
than even the actively star-forming Scd spiral galaxy tem- 



plate from CWW, and most are bluer than even the CWW 
Magellanic Irregular SED. The CWW Scd model would 
predict Vgoe ^ -f^ieo si2for2<z<3.5, redder than all 
but a very few of the NICIvIOS-selected galaxies in either 
the spectroscopic or photometric redshift based samples 
shown in Figure 1. Even for a sample selected from deep 
infrared images, it appears that nearly all HDF-N galaxies 
s-t 2 < z < 3.5 have specific SFRs (i.e., the instantaneous, 
ongoing rate of star formation relative to the integrated 
past star formation or total stellar mass) that are higher 
than those of present-day Hubble Sequence galaxies. Such 
rapid star formation is more characteristic of that in local 
starburst galaxies, which are not represented among the 
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CWW models. 

In Figure 2 we also compare the LBG photometry to six 
empirical starburst galaxy templates (Kinney et al. 1996) 
with varying dust extinction ranging from E{B — F) ~ 
0.0 — 0.7. The LBGs all fall within the envelope defined 
by these templates, which demonstrates that their UV- 
optical SEDs are broadly similar to those of local starburst 
analogs. The UV spectra are generally redder than those 
of the unreddened starbursts, and there is also a spectral 
inflection around the B aimer break region that indicates 
a strong contribution to the SED from longer-lived stars 
(A-type and later). 

We also fit all the Kinney et al. galaxy templates (five 
quiescent galaxy types - Elliptical, SO, Sa, Sb, Sc - six 
reddened starbursts, and the spectrum of the very blue, 
essentially unreddened starburst NGC 1705) to the SED of 
each LBG. The distribution of the best fitting templates is 
shown in Figure 3. The photometry for most of the LBGs, 
19/33 (^ 58%), are best fit by starburst templates with a 
moderate level of dust extinction, E{B — V) ^ 0.1 — 0.2. 
Of the other LBGs in the sample, 10/33 of the galaxies 
are well fit by essentially unreddened starburst models 
or the extreme blue template for NGC 1705. Only 4/33 
cases have best fit templates with large color excesses, 
E{B - V) > 0.25. None of the LBGs were weU fit by 
the quiescent galaxy templates. Based on these results, 
therefore, the typical LBG SED is broadly similar to local 
starburst galaxies with modest amounts of the dust ex- 
tinction — typical UV continuum suppressions at 1700 A 
on the order of, Anoo ~ 0.5 — 1.0 mag.^ These values 
are somewhat lower than those derived from other meth- 
ods (Sawicki & Yee 1998; Meurer et al. 1999; Steidel et al. 
1999, see also §5.1). 

One possible concern with the infrared LBG photom- 
etry might be the presence of a significant flux contri- 
bution from nebular emission lines, arising from heated 
and shocked gas associated with regions of star forma- 
tion (the most common being [O II] , [O III] , Ha+ [N II] , 
and H/3, which are labeled in Figure 2). Very strong 
emission lines can enhance the flux in a given band- 
pass and alter the galaxy colors. In several instances, 
strong nebular emission lines have been known to substan- 
tially affect broad-band infrared magnitudes and colors for 
high redshift radio galaxies (e.g., Eisenhardt & Dickinson 
1992, Eales et al. 1993), although generally radio galaxies 
have much stronger line emission than do ordinary, star 
forming galaxies. To investigate their contribution to the 
integrated fluxes, we consider the effect of an emission line 
with equivalent width W = Fune/Fx^ where F^ne is the in- 
tegrated flux in the emission line and F^ is the flux density 
of the continuum. Assuming a (roughly) flat continuum, 
the average flux density measured through a bandpass is 
approximately, F\ ~ -|- Fune AA~^, where AA is the 
width of of the bandpass. This can be reexpressed as 

:^ ~ 1 + + 



equivalent width Wq introduces a magnitude increase of 

I^o(l + ^) 



Am ~ —2.5 log 



1 + 



AA 



(2) 



For galaxies in the redshift range of interest here, the opti- 
cal nebular lines pass through the NICMOS and Kg fllters. 
The range of nebular emission lines observed in both the 
Kinney et al. local starburst templates and the Pettini 
et al. (2001) spectroscopic observations of 2; ~ 3 galaxies 
show that the equivalent width of [O III] can be as high as 
Wq 200 A. For the redshifts considered here, this corre- 
sponds to an increase in the observed broad-band fluxes 
of Am ~ 0.2 mag to the Hiqq- or i^g-bands. 

We investigated the emission-line contamination by ex- 
amining the residuals in the AFi/fx. distribution, where 
AFi = ,f{\i) — SiTj{\i) is the difference between the ob- 
served flux through a bandpass, i, and the "best-fitting" 
Kinney et al. template, j. In Figure 4, we show this distri- 
bution for the LBG sample. Each data point is coded to 
indicate if the bandpass contains one of the strong emis- 
sion fines, [O II], [O III], [N II], Ra, or H/3. Properly ac- 
counting for the photometric measurement uncertainties 
(which are generally larger for the infrared data), the scat- 
ter for bandpasses spanning one or more strong emission 
lines is equivalent to that for bandpasses devoid of emis- 
sion lines. To test this further, we refit the observed LBG 
photometry to the Kinney et al. templates with the major 
emission lines removed by interpolating across those wave- 
length intervals (i.e., [Oil], [O III], H/3, Ha, and [Nil]). 
The strongest emission lines in the local starbursts have 
similar equivalent widths to the largest values measured 
for LBGs by Pettini et al. (2001). The resulting distribu- 
tion of best-fitting spectra is essentially unchanged with 
respect to the previous results. Unless the nebular-line 
strengths for the HDF LBGs are significantly larger than 
those measured for other LBGs by Pettini et al., it seems 
unlikely that the emission lines have a significant effect on 
the spectral template fitting. 

4. COMPARISONS WITH STELLAR POPULATION 
SYNTHESIS MODELS 

Although we have seen that the spectral energy distribu- 
tions of LBGs are broadly similar to those of present day 
starburst galaxies, we must refer to models in order to 
characterize their stellar content and star formation histo- 
ries. Here we compare the LBG spectrophotometry to stel- 
lar population synthesis models to investigate constraints 
on the galaxy ages, star formation histories, dust extinc- 
tion, and stellar content. Because our photometry spans 
a broad range of rest-frame wavelengths (0.08-0.6 /xm at 
the mean redshift of the sample), we are sensitive to stel- 
lar population components with a wide range of ages and 
mass, and to different aspects of the galaxies' star forma- 
tion histories. Meurer et al. (1999) have demonstrated that 
the UV spectral slope (1000-4000 A) in nearby starburst 
galaxies correlates with the ratio of their far-IR to bolo- 
metric luminosities, and hence the degree of extinction for 

^ Kinney et aL derive the reddening in their starburst sample from the Balmer-Une decrement, the ratio of H|8-to-Ha: fluxes, which is 
a measure of the color excess in the nebular-emission regions, E{B — V)g. Calzetti, Kinney, & Storchi-Bergmann (1994) and Calzetti et 
al. (2000) have shown that there is a scaling relationship between this and the effective color excess of the stellar continuum [denoted as 
E{B — V)s]- Therefore, the color excesses of the Kinney et al. starburst templates correspond to a suppression of the UV continuum at 1700 A 
by A1700 ^ 4.6B(B - V)g mag. 



AA 



(1) 



where we substitute the rest frame equivalent width Wq = 
W{1 + z)^^. Therefore, an emission line with rest frame 
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Fig. 3. — Distribution of the starburst galaxy template types Kinney et al. (1996) that most nearly match galaxies in the LBG sam- 
ple. The columns correspond to each starburst template, as labeled. The first column refers to the very blue starburst galaxy, NGC 1705 
[E(B — y) ~ 0.0]. The next five columns are labeled by the range of reddening values that correspond to the averaged starburst templates. 



star-forming regions. The UV-to-optical flux ratio (e.g., 
the inflection around the Balmer/4000 A break) gauges the 
ratio of early-type (OB) to late-type (A and later) stars, 
which is a diagnostic of the past star-formation history. 

To generate model galaxy spectra, we used the newest 
version of the Bruzual & Chariot (1993) stellar popula- 
tion synthesis code (Bruzual & Chariot, in preparation, 
BC2000). The BC2000 models generate spectra of inte- 
grated stellar populations for a specified metallicity, ini- 
tial mass function (IMF), and star-formation history as 
a function of the population age. This technique repre- 
sents the galaxies' stellar content as the sum of a series 
of isochrones as a function of time for the specified star- 
formation history, using theoretical evolutionary tracks to 
predict effective temperatures and luminosities for stars 



as a function of mass and age. To compute the emergent 
spectrum, the BC2000 code uses stellar-spectral libraries, 
summing over the distribution of stars present at each time 
step. For a further description of the BC2000 models and 
references, see Bruzual (2000, and references therein). 

We will initially consider stellar populations formed by 
a single, continuous episode of star formation whose rate 
decays exponentially with a characteristic timescale, which 
we vary as a free parameter. The "age" of the galaxy is 
therefore defined as the time since the onset of the forma- 
tion of the stellar population that dominates the observed 
SED. The LBG star-formation histories might, in princi- 
ple, be significantly more complex than these monoton- 
ically evolving parameterizations, due to discrete events 
such as mergers, tidal interactions, shocks, or gas accre- 
tion. However, to model such stochastic processes would 
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Fig. 4. — Residuals between the observed LBG photometry and the best-fitting starburst galaxy templates. The data points show the 
fractional difference between the data and model versus the rest— frame wavelength, and error bars indicate the observed fiux uncertainty. 
Each data point is coded to indicate whether the bandpass spans a prominent emission line, as indicated in the figure. 
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require one to invoke an infinite number of discrete, possi- 
ble, star-formation histories. The choice of more complex 
star-formation histories would not greatly affect many 
of the conclusions here, which pertain to the nature of 
the stars that dominate the observed light in the ultra- 
violet and optical rest frame. Our definition of galaxy 
age depends to some degree on the form of the star 
formation history, because stars from earlier (hence older) 
star-formation episodes may remain undetected (see §5.2). 
The exponential, star formation histories provide a lower 
bound on the stellar masses of the LBGs, because they 
include possible models with very young ages, very active 
star formation, and hence the smallest values of M-jL. 
Thus, we should consider the stellar-mass estimates from 
these simple models strictly as lower limits. Below, in 
§5.2, we will consider the effects of a two-component 
star-formation history, in which the first component rep- 
resents the simple, monotonic, exponentially decreasing 
star-formation histories as above, and the second results 
from a stellar population formed in an instantaneous burst 
at 2; = 1000. Such a stellar component has a maxi- 
mal M-jL by design, and as such, we argue that the to- 
tal stellar-mass estimates from these models represent a 
conservative upper bound on the true stellar mass of the 
LBGs. Therefore, using these two sets of star-formation 
histories, we are able to bound the possible range of the 
total stellar masses of the LBG stellar populations in this 
sample. 

4.1. Model Parameters 

We generated a suite of population synthesis models 
spanning a wide range of the physical parameter space 
of age, star-formation timcscalc, mctallicity, IMF, and 
dust extinction. We allowed for several forms of the 
IMF (Salpeter 1955; Miller & Scab 1979; Scalo 1986), 
all of which adopt lower and upper mass cutoffs at 0.1 
and 100 M.q. Wc parameterized the star formation his- 
tory by adopting an instantaneous SFR of the form, 
^{t) ~ exp(— t/r), where t is the time since the on- 
set of star formation, and r is the characteristic star- 
formation, e-folding time. We considered values in the 
range r — 0.01 — 10 Gyr, in quasi logarithmic intervals, as 
well as a constant star-formation model (in effect, r « 00.) 
We tracked the time evolution of the model spectra in log- 
arithmic intervals for ages from 0.126 Myr to 20 Gyr. 

The BC2000 models offer the choice of two sets of stel- 
lar spectral libraries. One set, primarily based on Kurucz 
model atmospheres (see Lejeune, Cuisinier, &: Buser 1997), 
covers a wide range of metallicities, Z ~ 1/200 — ^ Zq. 
Another set uses empirical stellar spectra from the com- 
pilation of Pickles (1998), but is available only for solar 
metallicity. Wc compared results from the solar mctallic- 
ity model atmosphere libraries to those from the empirical 
spectral libraries, and found essentially no change in the 
results when fitting the broad band LBG photometry. 

The metallicities of LBGs are only loosely constrained 
from current data. Optical and infrared spectra of the 
gravitationally lensed z = 2.7 galaxy cB-58 (Pettini et al. 
2000; Teplitz et al. 2000a) suggest values ~ 1/4-1/3 Zq, 
and NIR spectra of several other LBGs suggest values con- 

^ Note that Calzetti et al. (2000) slightly updates the attenuation pai 
In particular, Calzetti et al. (2000) revises the ratio of total to select 



sistent with this (Teplitz et al. 2000b; Pettini et al. 2001). 
This is higher than measurements for damped Lyman a 
(DLA) systems at similar redshift, which generally have 
metallicities Z < 0.1 Zq (Pettini et al. 1997; Pettini 2000). 
One might reasonably expect higher metallicities for LBGs 
than for DLA systems at the same redshift because the for- 
mer are sites of active star formation, which may rapidly 
enrich the surrounding medium. Therefore, we have used 
the stellar atmosphere models so that we were able to in- 
vestigate the effects of varying metallicity on the model 
fitting results. However, it is important to remember that 
the theoretical stellar atmospheres at sub-solar metallic- 
ities are relatively untested by empirical data. Leitherer 
et al. (2001) have recently created a library of HST ultra- 
violet spectra of LMC and SMC O type stars, and find 
a fairly weak metallicity dependence for some UV line in- 
dices, but there is not yet much UV spectral information 
for later-type stars, nor data on the overall UV-to-optical 
SEDs of composite stellar populations at various metal- 
licities. In §5 below we will note various ways in which 
the results of our parameter fitting depend on the model 
metallicities. These should be regarded with some cau- 
tion, due to the current dearth of empirical calibration for 
the low-metallicity models. 

We included dust extinction using two different prescrip- 
tions: the Calzetti et al. (2000) attenuation curve inferred 
from local starburst galaxies, and the Pei (1992) parame- 
terization of the extinction curve for the Small Magellanic 
Cloud (SMC). The empirical Calzetti et al. starburst at- 
tenuation curve has a much "greyer" wavelength depen- 
dence in the near ultraviolet than does the SMC law. In 
part, this is believed to be due to geometric effects re- 
sulting from a mixed distribution of stars and dust in ex- 
tended galaxies. The classical stellar reddening curve for 
the SMC represents a combination of absorption and scat- 
tering terms. For a galaxy, photons that are scattered but 
not absorbed will eventually emerge. Pei (1992) provides 
separate parameterizations for scattering and absorption 
in the SMC extinction curve, and it might be appropri- 
ate to use only his absorption terms. Instead, we have 
retained the full SMC extinction curve with both scat- 
tering and absorption terms, precisely because it differs 
most substantially from the Calzetti relation. Our moti- 
vation is to diagnose the degree to which our results for 
LBG fitting depend on the extinction model. In this case, 
the steeper UV extinction of the classical SMC law pro- 
vides a useful test. In reality, however, we might expect 
dust extinction in LBGs to behave more nearly like the 
local starburst examples. Often, therefore, we will focus 
of our attention on results derived using the Calzetti et al. 
(2000) starburst extinction law.''' In the discussion that 
follows, we parameterize the effects of dust as the attenua- 
tion, A\ = E{B — V)k{X), in units of magnitudes, at some 
wavelength, A. At our UV reference wavelength, 1700 A, 
A;(1700A) = 9.65. 

4.2. Fitting the Synthesis Models to the LBG 
Spectrophotometry 

To fit the population synthesis models to the LBG 
photometric data, we converted the synthetic spectra to 

eters previously presented in Calzetti et al. (1994) and Calzetti (1997). 
extinction, Ry , to 4.05 ± 0.80, from 4.88 in Calzetti (1997). 
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bandpass averaged fluxes. The BC2000 code generates 
synthetic; spectra as specific luminosity density. £x , in units 
of solar luminosity per Angstrom per unit solar mass. We 
transformed these to luminosity per unit frequency and 
applied attenuation for dust extinction, 



L^{Xo,t,T, Axo,M) 



cm*{t) 



£A(Ao,t,r)10-0-*^^o, (3) 



where A = Xo{l + z) relates observed- and rest-frame 
wavelengths, and M is the total galaxy mass in stars. The 
BC2000 model spectra are normalized so that their total 
mass (gas + stars) is 1 Al©, and the stellar-mass fraction 
m*(t) is provided for each time step. From this we obtain 
the flux density, 

F^z, X,t,T,A,„M) = {l + z) ^-o(>^^^-A^) 

where Dl{z) is the luminosity distance (cosmology depen- 
dent), and tigm(2, A) is the wavelength dependent flux 
suppression due to the H I opacity of the IGM at wave- 
lengths shortward of Lyman a (Madau 1995). We mul- 
tiplied each spectrum with the (dimensionless) bandpass 
throughput function T^, which is the total throughput in 
the sense that it includes the response from the telescope, 
filter, and instrument optical assemblies. Thus, we obtain 
integrated synthetic, bandpass-averaged fluxes, 



{F^{z,t,T, Airoo,M)) 



jT,F4z,X,t,T,Ax,)du/i 



, (5) 



where we have parameterized the dust attenuation as mag- 
nitudes at (rest-frame) 1700 A, AnQQ. 

Once a model has been fit to the broad band galaxy 
photometry, its star-formation rate ^{t) = 5'oexp(— i/r) 
can be determined from the model normalization. This is 
the intrinsic rate, i.e., the physical rate that the stellar 
population model is forming stars at the observed age t. 
The dust correction has already been taken into account 
in the template fitting procedure. 

Due to the nature of the BC2000 models, only discrete 
values for the IMF and metallicity are available, and so we 
hold these parameters fixed at specific values. Although 
we tested a wide range of values, Z = 0.01 — 3.0 Zq, and 
several forms of the IMF, we concentrate our presentation 
here on models with Z — 0.2, 1.0 Zq and Salpeter or Scalo 
IMFs. These metallicity values are consistent with the lim- 
ited observational data on LBGs (see §4.1). As we will see 
below, the results of our analysis do not strongly favor any 
particular range of metallicity for HDF -N LBGs, although 
the choice of metallicity does sometimes affect the param- 
eters of the best-fitting models. To determine the best 
fitting models, we allow the variables t, r, AnoQ, and Ai 
to vary as free parameters. We generated a suite of ^ 10'^ 
synthetic models that span a wide range of the physical pa- 
rameter space {t ~ 0.001 — 20 Gyr, ^1700 — 0.0 — 6.7 mag, 
T = 0.01 - 10 Gyr and r ^ 1000 Gyr). 

We derived best-fit models for each LBG by comput- 
ing a statistic for each model. We defined ^oy each 
model set of synthetic data points, 

[f!^^ -{F^\z,t,T,A,roo,M)) 



(6) 



where the fi, are the observed photometric flux densities, 
and cr(/y) are the pure photometric uncertainties. The 
fTsys terms represent systematic errors that encompasses 
several uncertainties, which we discuss below. Although 
photometry is available through seven bandpasses for all 
objects (UsooB^^oVeoelsiiJiwHieaKs), we always exclude 
the F300W measurement because this point is well be- 
low the Lyman limit for all LBGs and suffers from severe 
IGM attenuation. The relation provided by Madau (1995) 
represents an average value of the IGM attenuation as a 
function of wavelength and redshift, but cannot account 
for variations in riGM(A, -j) along individual sight lines 
due to the stochastic distribution of absorbing clouds. For 
z > 2.5, the F450W bandpass extends significantly blue- 
ward of Lyman q, where variations in the Lyman a forest 
density may also affect the observed flux, and therefore 
we also excluded that band when fitting objects at those 
redshifts. By excluding these bandpasses, we prevent un- 
known statistical deviations in the IGM from influencing 
the fit. 

In principle, the expected value of Xo per degree of free- 
dom should be unity. In practice, we found that the min- 
imum values for the best fitting models were gener- 
ally larger than this. This implied that either our models 
incompletely describe the true SEDs of the LBGs, that 
the measurement uncertainties were larger than the sim- 
ple photometric errors, <j{fv), or that the parameter error 
distribution is not Gaussian (or some combination of these 
effects) . It is not surprising that the stellar population syn- 
thesis models are incomplete, because we are forced to use 
discrete parameter values for metallicity and IMF. We are 
also restricting the range of possible star-formation histo- 
ries to simple exponential models. Moreover, the popula- 
tion synthesis models might well have systematic inaccura- 
cies (see e.g., Bruzual & Chariot 1993, Chariot, Worthey, 
& Bressan 1996, Bruzual 2000), and there are almost cer- 
tainly uncertainties or inherent variations in the extinction 
properties (see Calzctti 1997, Calzetti et al. 2000). Finally, 
it is not unlikely that we have underestimated the flux un- 
certainties for the HST photometry. The values we have 
used here are based on the measured pixel to pixel im- 
age noise, but do not account for systematic errors in the 
measurements as implemented by the photometry software 
(e.g., issues of background determination, aperture posi- 
tioning, etc.), as well as additional systematic terms due 
to flat flelding errors, PSF mismatching, photometric zero- 
point uncertainties, etc. Many of these sources of system- 
atic error would introduce flux deviations (either errors or 
differences relative to a model) proportional to the source 
flux itself. Using the results of image simulations that were 
done to analyze the errors and detection efficiency in our 
HDF/NICMOS catalogs (see Dickinson et al. 2001), we 
included an additional error term (Jsys/ fv ~ 3.5% for the 
HST photometry. The i^s^band photometry was mea- 
sured using a different method (§2) whose uncertainties 
have been extensively tested and calibrated, and thus do 
not include this extra term. Even with this additional 
uncertainty, the reduced values were still larger than 
one, and we therefore included an additional 4% system- 
atic flux error for all bands in order to force the values of 
the best-flt x^ per degree of freedom to be of order unity. 
In general, including these additional uncertainty terms 
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does not substantially change the distribution of values 
in the model space or the results for the best-fitting model, 
but it does affect the mapping of the distribution to 
percentage confidence contours on the fitted parameters. 
Overall, it has the net effect of increasing the uncertainty 
on the model parameters. Therefore, we expect that the 
parameter constraints from our modeling are, if anything, 
conservative. 

4.3. Constraints on the Stellar Population Parameters 

In Table 3, we list the parameters for the stellar popu- 
lation synthesis models that best fit the observed photom- 
etry for each galaxy, i.e., the models with the minimum 
X^. Columns 1 and 2 hst the NICMOS ID number and 
redshift, respectively. For each LBG, there are four best 
models (rows) that show the results for two parameteri- 
zations of the IMF (column 4) and two metallicities (col- 
umn 3). The best-fit parameters (population age, star- 
formation e folding time, ^1700, and total stellar mass) 
are given in columns 5 — 8. All models listed in Table 3 as- 
sume the Calzetti et al. starburst extinction law, although 
we have also fit SMC extinction models, and will refer to 
these models in the following discussion. The instanta- 
neous SFRs associated with the best-fit model is shown in 
column 9. The minimum x^ corresponding to the best fit 
is given in column 10. 

The statistic can be used to construct confidence in- 
tervals around the best-fit parameters. With our suite 
of models, we explored the likelihood that other models 
produce equally good fits to the data. To derive this like- 
lihood given a best-fit model with Xqj we computed the 
probability that the correct model has a value x^ that is 
greater than the observed Xo- For each model, we obtained 
the difference in x^ between the best fit and every other 
model, Ax^ = (x^~Xo)- For each LBG, curves of constant 
Ax^ in the multi-dimensional parameter space translate 
directly to confidence distributions on the parameters. To 
make this transformation, we performed Monte Carlo anal- 
yses using the LBGs' fiux and uncertainty distributions 
(including systematic uncertainties). For each galaxy, we 
constructed 100 synthetic realizations of the data by mod- 
ulating the observed fiuxes with a random amount drawn 
from the (normal) error distribution (including both the 
photometric uncertainties and the systematic errors). We 
then refitted the models to the synthetic data sets and de- 
rived new best-fitting models. When compared to the 
original best-fitting model, the curve of constant Ax^ 
that encompasses some fraction of the synthetic "best- 
fits" corresponds directly to confidence percentages (equal 
to the contained fraction of synthetic best fitting models) 
on the parameters of best-fits to the true data. For 100 
Monte Carlo realizations, the uncertainty on the mapping 
of a curve of constant Ax^ to a confidence percentage is 
a{X)/X w 7.1% for parameter X. 

Figures 5(a-c)-7(a-c) present examples of the types of 
confidence intervals we derive on the fitted parameters, 
shown in two-dimensional projections of the parameter 
space. Each figure contains a series of panels showing the 
parameter dependences of population age versus attenu- 
ation, age versus total stellar mass, and age versus the 
characteristic star-formation time scale, for assumed val- 
ues for the metallicity and IMF. In each figure we also 



mark the lookback times to redshifts z = 5 and 1000. In 
the mass age confidence interval panels, we indicate the 
characteristic stellar mass for present-day galax- 

ies, 1.44 X 10" /i^ Mq, derived by Cole et al. (2001) 
from their measurement of the the local K band lumi- 
nosity function and the optical-IR color distribution of 
galaxies, assuming a Salpeter IMF. 

In most cases there is a range of models that all fit 
the data reasonably well. In Figure 8, we show exam- 
ples of model spectra, each from within the 68% confi- 
dence intervals, but with substantially differing parame- 
ters. With the available photometry, the models fit the 
observed SED very well in the rest-frame UV blueward of 
the 4000 A and Balmer breaks. The differences between 
the models become larger around and redward of the Kg- 
band [see especially the top panel, NIC 814 (HDF 3-93.0), 
where the redshift is such that the ilieo-band spans the 
Balmer break]. Thus, these models are not well differenti- 
ated. Deeper ii'-band data, and especially measurements 
at longer rest-frame wavelengths, e.g., from the Space In- 
frared Telescope Facility (SIRTF), should be quite effective 
in narrowing the range of possible stellar populations for 
galaxies at these redshifts. 

Some parameters are better constrained than others, 
and there are some degeneracies between the fitting pa- 
rameters for individual galaxies. For example, both dust 
and age redden the colors of the galaxies. While extinc- 
tion strongly reddens the UV portion of the spectrum, the 
Balmer break amplitude is relatively unaffected, and thus 
serves as a useful age indicator, constraining the relative 
number of early- (OB) to late-type stars (A and later). 
However, even with precise NICMOS photometry, we are 
able to derive only rather loose c;onstraints on age and red- 
dening, and the confidence distributions for those param- 
eters tend to be elongated along an axis of anticorrelation 
(see, e.g., Figures 5a and 7a). We note that this age- 
extinction degeneracy is considerably weaker when fitting 
with models using the SMC extinction relation. 

The best parameter constraints are the derived LBG 
stellar-mass estimates. Based on fitting the simple star- 
formation histories to the observed photometry for any 
given LBG, the 68% range of allowable stellar masses is 
constrained a{logM) < 0.5 dex, with median ~ 0.25 dex. 
In part, this is due to the high quality of the near-infrared 
data, which probes rest-frame optical wavelengths that 
are less sensitive to mass to-light ratio (Ai/L) variations 
due to population age or extinction. However, even in the 
rest frame S-band, A4/L can vary substantially for stellar 
populations with different ages. For a given galaxy, the 
range of possible, intrinsic values for M/L is restricted 
by the UV to optical SED fitting, which constrains the 
age and reddening for each galaxy. Although, we have 
seen, these parameter constraints are somewhat degener- 
ate, (e.g., for age versus extinction or metallicity), these 
degeneracies tend to cancel out in terms of their effect on 
the extrinsic A4/L values, i.e., the stellar mass over the 
emergent (after extinction) luminosity. E.g., for spectral 
models that fit the photometry for a particular galaxy, the 
younger models will have lower Ai / L, but require more ex- 
tinction, which suppresses the light and raises M./L. The 
reverse situation applies to older, less reddened model fits. 
In this way, the allowable range of M/L is limited, re- 
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Fig. 5a. — Spectral synthesis model fitting results for one galaxy, NIC 503 (HDF 2-903.0, z = 2.233), which is fairly typical of the spec- 
troscopic LBG sample. The star— formation history used here is parameterized as decaying exponentially with e-folding time r. As described 
in the text, these models provide a lower bound on the total stellar mass of the galaxy. In each panel, we show the equivalent 68% and 95% 
(orange and green shaded regions, respectively) confidence intervals on various quantities plotted against stellar population age, defined as 
the time since the onset of star formation: (a) 1700 A dust attenuation, (b) star-forming, e-folding timescale, and (c) total stellar mass. 
Each of the four panels corresponds to a different combination of the assumed metallicity and functional form of the IMF (as labeled). The 
best-fitting values (see Table 3) are indicated by the solid cross— hairs. The vertical dashed lines mark the elapsed time between the galaxy 
redshift and higher redshifts Zf = 5 and Zf = 1000. In (c), a line indicates the characteristic stellar mass of present— day L* galaxies (Cole et 
al. 2001). 



suiting in tighter constraints on the stellar mass for each 
galaxy (cf. Figures 5c, 6c, 7c). This method for estimating 
LBG stellar masses is essentially the same as that used by 
Brinchmann & Ellis (2000) in their analysis of field galax- 
ies at < z < 1. 

It is important to note that the stellar-masses deriva- 
tions are based on fitting the observed LBG photome- 
try to simple, constant or monotonically decreasing, star- 
formation histories. The results, therefore, are largely 
driven by the most luminous stellar populations (namely, 
the most recently formed stars). Thus, the models pre- 
dict best -fitting M./L ratios for the stars that dominate 
observed light, and because these models fit the light 
from the youngest, highest-mass stars, we must consider 
these M-jL values as lower bounds for the stellar com- 
ponent in these galaxies (neglecting changes to the form 
of the assumed IMF, see also §5.1.1). With these mod- 
els, we have relatively weak constraints on the presence of 
any additional, older stellar components (e.g., from prior 
episodes of star formation) , which might have substantially 
higher M-jL ratios due to a predominance of lower-mass 
stars. Therefore, we regard the stellar-mass estimates us- 
ing these simple star-formation histories strictly as min- 
imum values for the galaxies' total stellar mass. In §5.2 



below, we consider the effects on the derived stellar param- 
eters by adding a second, maximally old, stellar component 
to the simple models above. 

The stellar-mass estimates do systematically depend on 
the functional form of the IMF and to the metallicity of 
the population synthesis model. Generally, increasing the 
model metallicities shifts the stellar masses to higher val- 
ues. The size of this offset depends on the assumed form of 
the IMF. For a Scalo IMF, increasing the metallicity from 
Z = 0.2 Zq to 1.0 Zq shifts the derived stellar masses by 
< 0.1 — 0.2 dex. However, for a Salpeter IMF, changing 
the metallicity over the same range introduces a system- 
atic mass offset in most cases of ~ 0.3 — 0.5 dex. Further- 
more, changing from Salpeter IMF to Scalo IMF mod- 
els (for fixed metallicity) produces systematic stellar mass 
shifts of S{logM) = 0.0 to 0.6 dex (median ~ 0.2 dex) for 
Z = 0.2 Zq, and 0.0 to 0.3 dex (median ~ 0.1 dex) for 
Z = 1.0 Zq models. 

Two of the galaxies, NIC 989 (HDF 4-445.0, z ^ 2.500) 
and NIC 782 (HDF 4-316.0, z = 2.801), have very unusual 
distributions of parameter values. These galaxies have 
very narrow loci in the age-attenuation plots, favoring very 
young ages, extremely high extinction, and large stellar 
masses. These models appear to be somewhat unphysi- 
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Fig. 6a.— Same as Figure 5a but for NIC 1352 (HDF 4-497.0, z = 2.800), which has one of the flattest (bluest) UV spectral slopes 
the galaxies in the sample. 
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Fig. 6c.— 
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Fig. 7a. — Same as Figure 5a but for NIC 814 (HDF 3-93.0, z = 2.931), one of the reddest and apparently most massive LBGs in the 
sample. 
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cal, with instantaneous SFRs ~ 30,000 A4q yr^^. The 
best-fit model spectra deviate from the observed U^qq- 
and _B45o-band data by many standard deviations. This 
could suggest that the published spectroscopic redshifts 
for these objects are erroneous. For example, NIC 782 is 
apparently well detected at C/300, which seems unlikely for 
an object at z = 2.8. The photometric redshift analyses 
of Budavari et al. (2000) favor lower redshifts of z ~ 1.5, 
and z c:i 1.8 for NIC 989 and NIC 782, respectively (see 
also Fernandez-Soto et al. 2001). We cannot exclude the 
possibility that these two objects have unusual SEDs that 
are not well fit by the models we have used, but for lack of 
additional information we will not consider these objects 
further in the remaining discussion of the LEG sample. 

5. DISCUSSION 

The constraints on the LEG, stellar-population pa- 
rameters span a wide proportion of the available multi- 
dimensional parameter space. The intersection of the con- 



fidence intervals of all the LEGs in the sample generally 
shows that there is no region common to all the objects 
with any significance. There is some small overlap in the 
age-extinction confidence region (a compactification of the 
larger parameter space), which is common to all the LBGs, 
albeit with a very small likelihood (order of a few per- 
cent). Even so, the corresponding star-formation histories 
for these models vary for each object. We also note that 
genuine object-to-object variation is suggested from the 
large range of observed LEG SED colors, which is beyond 
what is expected from the photometric uncertainties and 
from simply varying only the dust content in the galaxies. 
Thus, this seems to suggest that there is genuine variation 
in the objects' stellar-population properties. However, it 
is difficult (if not impossible) to prove genuine variation ro- 
bustly amongst the LEGs in our sample given the current 
data. In this section, we use the general constraints on 
stellar population parameters for the LEG sample, which 
we use to investigate the assembly history of LEGs. 
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Fig. 9. — Composite probability distributions for the model parameters of dust attenuation versus population age for the whole LBG 
sample. The models used here assume the starburst extinction law. The contours represent the 68, and 95% confidence intervals averaged 
over the whole sample (orange and green shaded regions, respectively), and thus trace the distribution of allowable model parameters among 
the HDF— N LBGs. The points show the best— fit parameter values for each individual LBG, and are coded to indicate galaxies with z < 2.95 
and z > 2.95, i.e., below and above the redshift at which the Balmer break redshifts into the -ffieo filter. The cross-hairs show the (geometric) 
mean value for the distribution. Each panel shows results for a different combination of metallicity and IMF. 



5.1. Properties of the LBG Stellar Populations 

In Figures 9-13, we show combined confidence intervals 
for the entire sample of LBGs (excluding objects 989 and 
782 for reasons given above). These plots were constructed 
by averaging the likelihood distributions on the fitting pa- 
rameters for the individual galaxies, and thus show the 
favored (and disfavored) regions of the parameter space 
among the entire population. The best-fit values of the pa- 
rameters for individual galaxies are marked by points. Dif- 
ferent symbols code galaxies at redshifts below and above 
z = 2.95. For the higher redshift objects, the 4000 A 
break and 3646 A Balmer break move through the NIC- 
MOS F160W filter bandpass, leaving only the lower-S/N, 
i^s-band data to sample rest-frame optical light, and thus 
making the stellar-population constraints more uncertain. 
In the sections that follow, we discuss the implications of 
the model fitting for the stellar population properties of 
the LBGs. 

5.1.1. Initial Mass Function 

The best-fitting models using any of the IMFs (Salpeter, 
Miller & Scalo, and Scalo) are consistent with the data. 
The best-fit models using the Scalo and Miller-Scalo IMFs 
tend to favor younger ages and slightly lower attenuation 
values because the steeper slope at the high mass end of 
these IMFs results in spectra that are inherently redder. 



In a detailed spectral analysis of the gravitationally lensed 
galaxy cB-58 {z = 2.7), Pettini et al. (2000) show that 
the P Cygni profile of the C IV A 1549 line requires a sig- 
nificant population of stars with masses > 50 M.q, and 
is inconsistent with a truncated IMF or one with a very 
steep high-mass slope, such as the Miller-Scalo form. 

Varying the IMF lower mass cutoff would also change 
the derived stellar masses. The BC2000 models used here 
assume upper and lower mass cutoffs of A4u = 100 AAq 
and Ml = 0.1 Mq- For a Salpeter IMF, stars with masses 
less than 2 contribute typically less than 1% of the 
1700 A luminosity. For fixed 1700 A luminosity and fixed 
upper mass cutoff, the total stellar mass varies with the 
IMF lower mass cutoff and slope x, [N{M) oc M''^ dM] 
as 

^'°'°'\ \n{Mu/Mi) x^2. 

So, for example, for a Salpeter IMF (x = 2.35) with 
Ml = 1 Mq, the total stellar mass would be 39% the 
mass derived with Mi = 0.1 Mq. 

5.1.2. Metallicity 

Increasing the metallicity reddens the model SEDs. This 
can arise due to the combination of several effects. Higher 
metallicity stars have shorter main sequence lifetimes. 
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Fig. 10. — Composite probability distributions for the model parameters of star-forming, e-folding time versus population for the LBG 
sample, assuming the starburst extinction law. The contours and symbols are the same as in Figure 9. The dashed line indicates the 
one-to— one correspondence. 
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Fig. 11. — Composite probability distributions for the model parameters of stellar mass versus population for the LBG sample, assuming 
the starburst extinction law. The contours and symbols are the same as in Figure 9. The cross hairs show the (geometric) mean value for the 
distribution, combination. 
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Also, for a fixed age, higher metalhcity populations have 
lower effective teraperatures, and have stronger metallic 
absorption features. In general, the broad-band photom- 
etry provides little constraint on the LBG metallicities. 
However, the choice of metalhcity can have a significant 
impact on the other derived parameters. As can be seen 
from Figures 9 and 11, increasing the model metallicities 
from 0.2 Zq to 1.0 Zq has the effect of decreasing the 
required extinction while increasing the population age. 
This effect is seen across the entire range of metallicities 
we have explored, Z = 1/200 to 3 Zq. As described in §4.3 
above, there is some degeneracy between age and extinc- 
tion when fitting the models to the photometry. In this 
context, it appears that there is an interplay between the 
effects of metalhcity, age, and extinction on the galaxy 
colors. The choice of metalhcity changes the color of 
the model SED at fixed age, and the best-fitting mod- 
els required to fit the photometry of a particular galaxy 
then tend to slide along the age extinction anticorrelation. 
We caution again that the theoretical stellar model atmo- 
spheres used by the BC2000 models are relatively uncali- 
brated by observational data at low metalhcity, and there- 
fore any results that show a strong metalhcity dependence 
should be regarded with caution. We also note that this 
age-extinction degeneracy is strong for models using the 
starburst extinction law, but is much weaker when the 
SMC law is used. 

5.1.3. Star Formation History 

The past star-formation histories for individual galaxies 
are only loosely constrained by the current data set, (see 
the age-r panels in Figures 5b-7b, and also Figure 8). 
Generally, a large region of the parameter space is con- 
sistent with the data, especially for t/r < 1. where the 
SFR is roughly constant and most values provide good 
fits. In Figure 10, we show the mean age-T confidence 
region for the LBG population. In most cases, models 
with t ^ T, where the current SFR is much smaller than 
the past average, arc strongly disfavored, although a few 
galaxies can be fit by "post-starburst" models with ages 
^100 Myr (see §5.1.6). In other words, the LBGs do not 
appear to be UV-bright objects due to a relatively small, 
residual tail of on-going star formation in an older galaxy, 
but rather arc seen when they arc still in the process of 
forming a significant fraction of their stellar content. This 
result may, in part, be due to the choice of exponential, 
monotonically decreasing or constant star formation his- 
tories; we consider other possible star-formation histories 
in §5.2 below. Although for most LBGs the star-formation 
timescales are only poorly constrained, we do note that for 
all combinations of metalhcity and IMF there are a num- 
ber of galaxies that are best fit by short timescale "bursts" 
(see Figure 10). 

5.1.4. Dust Extinction 

For models using the Calzetti et al. (2000) starburst ex- 
tinction law, there is some degeneracy in the attenuation- 
age parameter confidence regions (see discussion below). 
For individual LBGs, the 68% confidence region typically 
spans AAnoo ~ 1.5 — 2 mag, which corresponds to a range 
of allowable flux corrections that vary by a factor ~ 6 at 
1700 A. For the population as a whole (see Figure 9), the 



distribution of best-fit values spans -Ax7oo — — 4 mcig 
(~ - 3 mag) for Z = 0.2 Zq (1.0 Zq), with only a mild 
dependence on the chosen IMF, but with notable sensitiv- 
ity to metalhcity. For low metalhcity models, a number 
of galaxies arc best fit with young ages and heavy red- 
dening, but this is not seen when solar metalhcity models 
are used. The probability weighted mean UV attenuation 
for the LBG sample is {Anoo) = 1-6 mag (1.2 mag) for 
Z = 0.2 Zq (1.0 Zq). The average and range of Anoo 
values derived here are similar to those derived by Steidel 
et al. (1999) and Adelberger & Steidel (2000) based on 
optical photometry for a large data set of several hundred 
LBGs, using the same dust attenuation law but simpler 
stellar population assumptions. However, the luminosity- 
weighted average UV flux correction (i.e.. the ratio of in- 
trinsic to emitted luminosity summed over the whole sam- 
ple) is considerably larger, a factor of ~ 13.8 (4.6) for 
Z = 0.2 Zq (1.0 Zq) and a Salpeter IMF (and similar 
for the Scalo IMF). This reflects the tendency for LBGs 
with the largest fitted dust corrections to have the highest 
intrinsic UV luminosities and SFRs, a point emphasized 
by Meurer et al. (1999) and Adelberger & Steidel (2000). 
Those authors argue that this is a real trend among LBGs, 
and not an artifact of the extinction correction process. 

If the SMC extinction law is adopted instead, the anti- 
correlation between age and attenuation largely disappears 
(Figure 12), and the tail of young and high-extinction 
models found with the "grey" starburst dust attenuation 
is absent. The allowed model extinction at 1700 A is 
< 1.5 mag regardless of IMF and metalhcity, and any cor- 
relation with population age is weak at best. Most galaxies 
are well fit by stellar population ages > 40 Myr. 

Regardless of metalhcity, IMF, and dust model assump- 
tions, models with very young ages and low dust content 
are strongly disfavored. Very young, unreddened star- 
forming regions, particularly those with low metalhcity, 
should have UV spectra that rise (in units) toward 
shorter wavelengths. In practice, this is never observed 
among the LBGs: even the youngest galaxies apparently 
contain significant amounts of dust at these redshifts. This 
may suggest that the LBGs have experienced metal enrich- 
ment from a previous stellar generation, and therefore the 
LBG population probably are not "first-generation" ob- 
jects (unless the spectra of young, low metalhcity, dust- 
free objects are sufficiently different from the models used 
here). Such first- generation objects must be quite rare 
at these redshifts, or must pollute themselves with dust 
on very short timescales. Alternatively, perhaps the SFRs 
(and hence luminosities) of such objects are too low to 
allow them to be detected in the the HDF images. 

The extinction correction factors that we derive from 
the model fitting differ somewhat from those obtained 
by other authors using various techniques. Modelling 
many of the same HDF-N galaxies using ground-based 
NIR photometry, Sawicki & Yee (1998) find SED fits 
with somewhat higher median UV-flux-suppression fac- 
tors, {Anoo) ~ 2.7 mag for Z = 0.2 Zq models (ac- 
counting for small differences in the forms of the star- 
burst attenuation law). The two analyses are in closer 
agreement at solar metalhcity, where Sawicki & Yee find 
(^i70o) — 1-4 — 1.9 mag, depending on the assumed star- 
formation history. The differences between these analyses 
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Fig. 14. — Comparison of the fitted SED model extinction (using the starburst attenuation law), in magnitudes at 1700 A, with the 
attenuation inferred from the U-spectral index, alone {fx oc A'') using the relation from Meurer et al. (1999). The /3 values measured 
from the WFPC2 photometry are indicated along the top axis. The points and error bars in the ordinate show the best-fit values and 68% 
confidence intervals on dust attenuation from the population synthesis models fit to the optical— to-infrared photometry. 



may be due to several factors, including the improved sen- 
sitivity of our NICMOS photometry, more accurate match- 
ing of optical to infrared photometry, allowance for a wider 
range of possible star-formation histories (Sawicki & Yee 
considered only 5-function bursts and continuous star- 
formation models), and the use of a larger galaxy sam- 
ple. Meurer et al. (1997) also found large UV-attenuation 
values (~ 15 — 20) by fitting the the UV slopes of a sam- 
ple of HDF-N LBGs to the starburst extinction law from 
Calzetti et al. (1994). However, using a revised method 
of a locally calibrated UV slope to far-infrared (FIR) dis- 
tribution, Meurer et al. (1999) derive more modest mean 
extinction of factors ~ 5. As noted above, the attenua- 
tions we derive are similar to those found by Adelberger & 
Steidel (2000) using UV rest-frame data for a larger LEG 
sample. Overall, the fact that the derived UV-flux at- 
tenuation depends so strongly on parameters on which we 
have little reliable information, such as stellar-population 
metallicity, leads us to believe that one must regard the 
quantitative results of such analyses with considerable cau- 
tion, but the necessity for significant dust attenuation in 
many or most LBGs seems secure. 



In local, UV-bright starburst galaxies, the UV-spectral 
slope P (/a oc A^) correlates with the flux emitted at 
far infrared wavelengths (8-1000 /im), demonstrating a 
close relation between UV extinction and dust reprocess- 
ing. Meurer et al. (1999) calibrate this relationship as 
^1700 = 4.30-1-1.93/3 (where we have adjusted their param- 
eters slightly to account for a small shift in the reference 
UV wavelength). In Figure 14, we compare the atten- 
uation derived from the UV-through-optical rest-frame 
SED model fitting to that computed from the UV slope 
(3 alone, using the prescription of Meurer et al. (1999) to 
convert Vgoe — -^814 colors to /?. For a given set of IMF 
and metallicity assumptions, there is considerable scat- 
ter (~ 1 mag) in the comparison of attenuation values. 
This is appreciably larger (by approximately a factor of 
two) than uncertainty that results from random and sys- 
tematic errors in deriving (3 from UV-to-FIR flux ratios 
(see Meurer et al. 1999). The uncertainty in ^1700 can 
produce substantially differing UV-flux correction factors 
(by an order of magnitude in flux), which directly affects 
the inferred SFRs (see discussion below). In some cases, 
galaxies with f3 values that predict ^1700 — 0.5 — 2.5 mag 
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can be fit by models requiring negligible extinction. For 
these objects, the full SED model fitting would suggest 
that the UV-spectral slope is due to an aging stellar pop- 
ulation rather than dust. There are also systematic off- 
sets in the derived attenuations, depending mainly on the 
choice of metallicity for the population synthesis models. 
The lower metallicity models are bluer, and thus require 
more reddening to match the observed colors. On aver- 
age, the best-fit values of A^oo for the Z = 0.2 Zq syn- 
thesis models are greater than those predicted from the 
UV slope alone by S{Anoo) ~ 0.35 mag (0.20 mag) for 
Salpeter (Scalo) IMF, while the solar metallicity models 
nearly always result in attenuation values lower than the 
f3 predictions by ~ 0.6 mag for both Salpeter and Scalo 
IMFs. Therefore, we emphasize that while it may be jus- 
tifiable to apply an average flux correction to a galaxy 
ensemble, there are very likely large uncertainties, both 
random and systematic, when making extinction correc- 
tions to individual objects. 

5.1.5. Stellar Population Age 

As we have already noted, there is some degeneracy be- 
tween age and extinction when fitting models to the broad- 
band LBG photometry (see, e.g., §4.3 and Figures 5a and 
7a). A similar trend is seen for the composite distribu- 
tion in Figure 9. It is not clear how much of this anti- 
correlation can be explained by the combined effect of the 
parameter degeneracies for individual galaxies, or whether 
there is some real tendency for younger galaxies to have 
more extinction. This might not be unexpected, because 
young star formation within our own Galaxy tends to be 
highly dust obscured, and at later times, stars may blow 
away dust or migrate from their dusty birthplaces. On 
the other hand, one might expect very young galaxies to 
be relatively free of dust and metals, unless these were pro- 
duced during previous episodes of star formation, before 
the birth of the current generation of stars. 

As can be seen from the composite probability plots 
(Figure 9), the best-fit solar metallicity models span an 
age range of 30 Myr ^t<l Gyr. The upper range of pop- 
ulation ages is similar for both low {Z = 0.2 Zq) and high 
(solar) metallicity models, but the low metallicity models 
exhibit a tail of solutions with very young ages (down to 
~ 4 Myr) and high extinction, A1700 ~ 3 to 4 mag with the 
Calzetti ct al. dust model. As discussed above in §5.1.2, 
the low metallicity models are intrinsically bluer, and thus 
require more extinction to fit the observed colors at a fixed 
age. It is notable that the fitting solutions with very young 
ages and high extinction appear only when using the star- 
burst attenuation law, which has a rather "grey" wave- 
length dependence. Fits with the steeper SMC dust law 
avoid this young and heavily reddened region of parameter 
space. 

Although the population fitting with low metallicity 
models allows very young ages (< 10^ yr), these timescales 
may be somewhat unphysical based on the kinematics and 
linear size scales observed for LBGs. Nebular-emission- 
line widths of (J ~ 60 — 120 km s^^ have been measured 
for several LBGs (Pettini et al. 1998, 2001; Tephtz et al. 
2000a; Moorwood et al. 2000). The objects in our sam- 
ple typically have physical diameters of 3 — 10 kpc. At- 
tributing the emission-line widths to virialized motions. 



the corresponding dynamical timescales for these galax- 
ies are t^yn ~ 10'^'^~^ yr. This may arguably set a lower 
bound on the time required for star formation to propagate 
across the galaxies, and thus on the age of their stellar pop- 
ulations. Therefore, these solutions are possibly disfavored 
even if formally allowed by the SED fits, unless one invokes 
an unusual geometry in which the star-formation regions 
and nebular emission lines do not originate in identical re- 
gions. In §5.1.6 below, we consider constraints on these 
solutions from limits on sub -mm emission in the HDF N. 

The geometric mean age for the composite distribution 
is ~ 70 Myr (120 Myr) for models with metalHcity 0.2 Zq 
(1.0 Zq) and Salpeter IMF. For Scalo IMF models, the 
distributions have slightly lower mean ages, ~ 40 Myr 
(70 Myr) for 0.2 Zq (1.0 Zq) metalhcity. For the solar 
metallicity, Salpeter IMF model fits, three galaxies have 
best-fitting model ages t > 1 Gyr, resulting from their 
moderately red -ffieo — K colors and model fits implying 
little or no dust extinction. All three galaxies, however, 
are at z > 3.2, where only the Ks data reaches rest frame 
A > 4000 A, and for one object (NIC 284, by far the 
faintest object in the sample), the photometry provides 
only a weak upper limit to the Ks flux. For this reason, 
the possible older ages for these galaxies must be regarded 
with caution. 

Sawicki & Yee (1998) derived upper limits oi t < 
200 Myr on the ages of the dominant stellar populations in 
most HDF LBGs, and concluded that LBGs are dominated 
by very young populations [t ~ 25 Myr). Their models 
with (5-function star-formation histories have UV lumi- 
nosities that decline rapidly, mandating young ages for a 
UV-bright (by definition) LBG population. It is there- 
fore more useful to compare their constant star formation 
models to our results. As can been seen from the indi- 
vidual galaxy confidence distributions and the composite 
age-attcmiation plots, many of the LBGs in our sample 
have allowable models with ages on the order of t ~ 1 Gyr. 
However, it is also true that for the majority of galaxies, 
our best-fit ages for low metallicity, Salpeter IMF mod- 
els are also < 200 Myr, with a weighted mean value of 
~ 70 Myr, only ~ 2x the median age found by Sawicki 
& Yee for z < 3 galaxies assuming the same metallicity 
and constant star formation. Adopting more metal rich 
models or a steeper UV-extinction law shifts the age dis- 
tribution to larger values, with median or weighted mean 
values in the range of 100 — 300 Myr and a substantial tail 
extending to ~ 1 Gyr or more. Although we generally find 
somewhat older ages than do Sawicki & Yee, our results 
do not exclude very young ages of several tens of Myr for 
many objects, at least not for fits using the low metallicity 
models. 

The ages derived from these best-fit stellar population 
models are generally much younger than the age of the uni- 
verse at the LBG redshifts. In Figure 15, we plot the 68% 
confidence interval on the model ages for each LBG ver- 
sus the observed galaxy redshift. We also overplot curves 
showing the lookback times to higher redshifts. The re- 
sults are somewhat dependent on the particular choices 
for the IMF and metallicity, and derived ages only apply 
to the youngest stellar populations under the chosen star- 
formation histories, which do not necessarily apply to the 
galaxies as a whole (i.e., the star-formation history for the 
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Redshift Redshift 

Fig. 15. — The distribution of LBG stellar population ages as a function of redshift. The data points and error bars show the best-fit ages 
and their 68% confidence interval from the stellar population synthesis models. The dashed curves show the lookback time from the observed 
galaxy redshift to constant "formation" redshifts, zj = 2.0, 2.5, 3.0, 3.5, 4.0, 5.0, 10.0, and 1000 (not all curves are labeled). 



entire galaxy may indeed be more complex than the simple 
monotonically decaying models used here). The derived 
ages from these models generally favor a scenario where 
most of the galaxies have formed their observed stellar 
populations (i.e., those stars that dominate the rest-frame 
UV-to-optical SED) within a rather small redshift interval 
prior to the redshift at which they are observed, at least for 
those objects at z < 3.0, where the model age constraints 
are best. If this is the case, then these LBGs assembled 
(or, are assembling) their stellar masses relatively rapidly. 
The LBG sample is devoid of old (i > 1 Gyr), quiescent 
galaxies at 2 < z < 3, which would be expected to popu- 
late the upper-left region of the panels in Figure 15. This 
might be expected, because the galaxies are selected to 
be UV-bright, but as we have previously noted, there are 
few (if any) candidates for old, quiescent galaxies among 
the NICMOS infrared-selected HDF-N catalog with pho- 
tometric redshifts 2.0 < Zph < 3.5 (see Figure 1). I.e., 
there is little evidence for a population of "fading," post- 
burst remnants from evolved galaxies formed at z > 3. 
However, the presence of older stellar components within 
the star-forming LBGs, cannot be ruled out. As we will 
see in §5.2, such older stars could, in principle, provide 
a significant fraction of the galaxies' total stellar masses. 
Under such a scenario, the fitted ages for the young stellar 
components might not represent the ages of the dominant 
(by mass) stellar populations. 



As with extinction, the strong dependence of the de- 
rived age distribution on relatively uncertain parameters 
such as metallicity, the extinction law, and the assumed 
prior star-formation histories means that one must inter- 
pret these results with caution. The question of LBG ages 
cannot be completely resolved with the available broad- 
band photometry, and would benefit from a more detailed 
knowledge and independent data on the chemical enrich- 
ment history and dust content of these galaxies, as well as 
from measurements of starlight at still longer rest-frame 
wavelengths. 

5.1.6. Star Formation Rate 

In Table 3, we list the instantaneous star-formation 
rates for the best-fit models for each LBG. One interesting 
result is the wide range of SFRs [10"^ < \E'/(Al0 yr~^) < 
10^] derived for galaxies in the sample, and indeed the 
wide range even among acceptable models for individual 
galaxies. The highest derived SFRs mostly correspond to 
models with young ages (mostly < 15 Myr) and small t/T. 
Most of these are the young, low metallicity, high extinc- 
tion models discussed above. For ages < 10^'^ yr and small 
t/r, the UV luminosity emitted per unit SFR is lower than 
its asymptotic value at later times because the OB stars 
that produce most of the UV light are still building up on 
the main sequence. At later times, their birth and death 
rates equilibrate, and the UV luminosity per unit SFR sta- 
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Fig. 16. — Predicted 850 fim flux densities from the LBGs based on the stellar population and extinction model fitting. The sub— mm fluxes 
were predicted using the FIR-UV flux density correlation in local starburst galaxies (Meurer et al. 1999; Adelberger & Steidel 2000), with 
the dust attenuations derived from the best— fitting models. The four panels correspond to various combinations of IMF and metallicity, as 
labeled. The yellow shaded region in each panel indicates sensitivity range of the SCUBA 850 fim survey of the HDF-N, /y{850^m) > 2 mjy 
(Hughes et al. 1998). For the 0.2 Zq models, one might expect that a few of the youngest, dustiest, and most actively star-forming galaxies 
could have been detected with SCUBA. 



bilizes. This, combined with the large extinction imphed 
for these models, results in model SFRs that are substan- 
tially larger than those which would be derived from the 
UV-luminosity calibrations generally used in the literature 
(e.g., Madau, Pozzetti & Dickinson 1998). The large SFRs 
are required in order to produce the inferred stellar mass 
on such short timescales. Conversely, a few galaxies can 
be fit with negligible amounts of ongoing star formation. 
These are "post-starburst" models with t ^ t and ages 
90-150 Myr, where the UV light primarily originates from 
B stars that have not yet left the main sequence. 

For the young, low metallicity models with high extinc- 
tion, the derived instantaneous SFRs are extremely high, 
(SFR > 1000 Mq yr~i). With this combination of high 
SFRs and extinction, one would expect that most of the 
energy from star formation should be absorbed by dust 
and reradiated in the far-infrared, perhaps with signifi- 
cant emission at mid-infrared and radio wavelengths as 
well. To date, however, there have been no robust detec- 
tions of LBGs in the HDF-N at mid-IR, far-IR or sub-mm 
wavelengths (Goldschmidt et al. 1997; Hughes et al. 1998; 
Downes et al. 1999, although there are two possible radio 



detections, Richards 1998). Therefore we can potentially 
use these surveys to constrain the intrinsic SFRs of these 
LBGs. 

In particular, we considered the SCUBA 850 /im ob- 
servations of the HDF-N (Hughes et al. 1998), because 
sub-mm observations are, in principle, quite sensitive to 
star formation in galaxies from z ~ 1 — 10 due to the 
strong negative fc-correction. Using the relations given by 
Adelberger & Steidel (2000), we predicted LBG flux den- 
sities for SCUBA at 850 /xm from the observed LBG (rest- 
frame) UV luminosity and the dust attenuations from the 
best fit models. Figure 16 shows the predicted 850 /xm 
SCUBA flux density as a function of the LBG stellar pop- 
ulation age. For the solar metallicity models, all of the 
LBGs have predicted 850 /xm flux densities that are below 
(or very close to) the 2 mJy hmit of the HDF-N SCUBA 
observations. For the Z = 0.2 Z© models, however, a few 
of the galaxies would have predicted 850 fim flux densities 
> 2 mJy. In most cases, the predicted sub-mm emission is 
only slightly above the SCUBA detection threshold. Given 
the range in the allowable Anoo values for each galaxy, the 
lack of SCUBA detections may not conclusively rule out 
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any particular model.® Moreover, recent mm and sub-mm 
observations of gravitationally Icnscd high redshift galaxies 
(van der Werf et al. 2000; Baker et al. 2001) apparently find 
less far-IR emission than predicted from the far-IR/UV 
correlations derived at z = 0. Conversely, Meurer ct al. 
(2000) find that the UV FIR correlation for local star- 
bursts substantially underpredicts the far-infrared emis- 
sion from nearby ultraluminous infrared galaxies. At this 
time, then, it seems difficult to set additional constraints 
on the SFRs derived from our SED model fitting. Ulti- 
mately, if the SFR can be constrained from observations 
in the FIR, then this could constrain the viable population 
synthesis models for each LEG. 

5.1.7. Stellar Mass 

Figures 11 and 13 show the distribution of best-fitting 
stellar-mass estimates for the spectroscopic LEG sample, 
under the assumption of the Calzctti and SMC extinction 
laws, respectively. Recall that due to the simple nature of 
the star-formation histories of the models, we treat these 
stellar mass derivations strictly as lower limits; we will 
consider upper limits to the stellar masses below (§5.2). 
It is also important to note that this is not a complete 
sample of galaxies, and certainly not a mass-limited one. 
The distribution of masses reflects this fact, e.g., in terms 
of the lower mass limit for the sample. 

The (geometric) mean stellar mass of the composite dis- 
tribution is ~ 3(6) X lO'^ Mq for the Z = 0.2 Zq (1.0 Zq) 
metallicity and Salpeter IMF models, assuming the star- 
burst extinction relation. Using a Scalo IMF shifts this 
value to =i 7 X 10^ AAq (insensitive to metallicity). The 
68% confidence range of our composite model mass dis- 
tribution spans the interval M. ^ 10^^^^ A^o, similar to 
the range of best-fit values for individual LEGs. The in- 
ferred stellar mass would be quite sensitive to the lower 
mass cutoff of the IMF, a parameter that we do not vary 
here (see §5.1.1). Thus, our conclusions are based on the 
assumption that this cutoff does not change from galaxy 
to galaxy. There is a weak correlation between the stellar 
mass and the population age for the Salpeter IMF models 
(no correlation is apparent for either the Scalo or Miller 
& Scalo IMF models). The distribution of stellar masses 
for the LEG sample is nearly independent of the assumed 
dust-extinction law (cf. Figure 11 and Figure 13). We 
attribute this primarily to the high-quality near-infrared 
data, which are sensitive to the later-type stars that con- 
stitute the bulk of the stellar mass, and where the effects of 
dust extinction are reduced. The LEG stellar masses are 
nearly all smaller than those of present day, L* galaxies 
(1.4 X 10"/ito Mq for a Salpeter IMF, Cole et al. 2001), 
but are similar to that of the Galactic bulge, ^ 10-'^° M.q 
(Rich 1999). Sawicki & Yee (1998) estimated similar me- 
dian masses, 3 — 8 x 10^ A^©, for Z = 0.2 Zq, Salpeter 
IMF models (transforming to our assumed cosmology) . In 
a new, infrared photometric study of a sample of more lu- 
minous (mostly L > L^y) LEGs, Shapley et al. (2001) 
derive median stellar masses (for Z = Zq, Salpeter IMF 
models) that are ~ 4x larger than those of our fainter 
HDF objects. 

^ The z = 2.93 galaxy NIC 522, or HDF 1-54.0, is the only LBG with y. 
solutions, with /^(850;um) ~ 3 - 10 mjy (~ 1 - 3 mjy) for Z = 0.2 % 
UV spectral slope, and we note in passing that this object lies only : 
[/^(850//m) = 3.0 mJy] of Hughes et al. (1998). 



Somerville, Primack, & Faber (2001) have considered 
various models for the nature and evolution of high red- 
shift galaxies. In particular, they present distributions of 
stellar ages and masses for "mock-HDF" catalogs derived 
using their "coUisional starburst" model, where high red- 
shift star formation is driven primarily by mergers and 
tidal interactions. For bright (Vgoe < 25.5) galaxies in the 
redshift interval 2.0 < z < 3.5, Somerville ct al. predict 
a broad range of stellar masses, 10^ < M/Mq < 10^^, 
with a median value of M ~ S x 10^ Mq (for h = 0.7). 
This is somewhat larger than the mean value (for Salpeter 
IMF models, as adopted by Somerville et al.) from our 
small HDF-N sample, but the comparison is not precise 
because our sample is not complete and magnitude lim- 
ited at Vgoe < 25.5. The Somerville ct al. models include 
a small tail of objects with masses greater than those of 
present-day L* galaxies, ^A > 10^^ Mq, which are not 
seen in the HDF N sample. However, this may simply 
be due to small number statistics of the HDF-N itself. 
More notably, the coUisional starburst model also predicts 
a significant number of bright (Veoe < 25.5), low-mass 
(10^ < M/Mq < 10^) galaxies, which are also not found 
in our LBG sample. These arc presumably small galaxies 
that have recently experienced a massive episode of star 
formation, which pushes them above the optical magni- 
tude limit. They could also be young, dust-free objects, 
which are also not found in our SED analysis. A compar- 
ison with HDF N LBGs selected from a Vgoe based cat- 
alog shows that none with Vgoe < 25.5 are missed in the 
NICMOS-selected sample. Therefore, given that the col- 
ors and fluxes of our spectroscopically selected LEGs are 
representative of all galaxies with photometric redshifts in 
this interval in the HDF N (see Figure 1), the absence of 
optically bright but low-mass LEGs may point to a dis- 
crepancy with the coUisional starburst models. 

We find a fairly good correlation between the derived 
stellar masses and the rest frame UV luminosity in these 
galax;ies (see Figures 17 and 18). This result is somewhat 
surprising, because the UV light is generally expected to 
trace the instantaneous SFR but not necessarily the stellar 
mass. The scatter of the best fit values in this correlation 
is (7(\ogM/ Ltjv) ~ 0.3 dex, regardless of IMF and metal- 
licity. This correlation could be misleading if the full range 
of HDF-N galaxies with 2.0 < ^ph < 3.5 do not follow this 
trend. This would require that the low-mass end of the 
photometric redshift galaxy sample contain a large popu- 
lation of old or highly extincted galaxies, which would pro- 
duce larger Veoe— -f^ieo colors at fainter magnitudes. There 
is no indication of any such population; indeed, if anything 
the LEG color-magnitude diagram. Figure 1, shows that 
most of the redder LBGs are brighter. In the next section, 
we will consider limits on "hidden" stellar mass from pre- 
vious generations of star formation, and it may be that 
the apparent L^/y-mass correlation would be weakened if 
such mass were indeed present. 

The stellar masses of the brighter LBGs, with UV lumi- 
nosities near the L* value from Steidel et al. (1999), are 
of order 10^° Mq. It is interesting to compare this with 
kinematically estimated masses based on measurements of 

icdicted 850 /lim flux density well above the SCUBA limit for all model 
(1.0 Zq). This object has a clearly disturbed morphology and a steep 
; 7" from the reported centroid position of the SCUBA source, 850.3 
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Fig. 17. — LBG stellar masses versus observed, rest-frame UV luminosity. The different panels show results from exponential star-formation 
history models with various metallicity and IMF assumptions, and all assuming the Calzetti et al. (2000) starburst extinction law. The points 
mark the best-fit masses, with error bars indicating the 68% confidence interval. The UV luminosity is given as an absolute magnitude at 
1700 A on the AB scale. The characteristic UV luminosity L*i7oo for LBGs at z « 3 from the luminosity function of Steidel et al. (1999) 
corresponds to M*^qq = —21.0, marked by a vertical dashed line, while the stellar mass of a present-day L* galaxy is indicated by the 
horizontal dashed line. In each panel, the diagonal dashed line corresponds to a constant mass-to-light ratio M/L\jy. The solid lines show 
a least squares fit of a straight line to the the data. 



nebular-emission-line widths for LBGs (see §5.1.5). A 
naive application of the virial theorem combining these 
line widths and the typical half light radii for LBGs sug- 
gests masses that are also ^ 10^° ■M-q- Thus, the stellar 
masses that we infer for ^ L* LBGs can account for most 
or all of the kinematically estimated mass. In fact these 
stellar-mass estimates are, if anything, lower limits, mea- 
suring only the mass of the "young," actively star form- 
ing component of the galaxy. Furthermore, as described 
in §5.2 below, the total stellar mass could be larger if a 
maximal old stellar population were present. As Pettini 
et al. (1998, 2001) have emphasized, the nebular emission 
lines probably originate in the UV bright, central regions 
of the galaxies, and thus may not sample the full gravita- 
tional potential well. It is therefore unclear whether these 
virial estimates really measure the total mass of the LBG 
dark matter halos. Unless the ratio of dark to luminous 
matter in LBGs is much smaller than that in present day 
galaxies, it appears that the simple virial mass calcula- 
tions based on the nebular line widths must significantly 
underestimate the total masses of the dark-matter halos 



in LBGs. 

5.2. Constraints on old stellar populations at 2 < z < 3.5 

The masses computed from the best fitting SEDs repre- 
sent the integrated stellar mass for the stellar population 
that dominates the rest frame UV-to-optical light of the 
galaxy. Strictly speaking, these estimates constitute only 
a lower limit to the total integrated stellar mass. We may 
speculate about the existence of older stars in these galax- 
ies that do not substantially contribute to the observed 
light, but which might contain a significant fraction of the 
total stellar mass. Observations of star clusters in the 
Galactic bulge (e.g., Mc William & Rich 1994, Ortolani et 
al. 1995) have shown that ages of their stellar populations 
are consistent with the age of the Universe, implying very 
early formation epochs (see also Renzini 1999). In princi- 
ple, a significant mass fraction of old stars could be present 
in LBGs, but invisible at UV to optical wavelengths, hid- 
den under the glare of younger stars. 

In §2 we noted that nearly all infrared-selected HDF- 
N galaxies with spectroscopic or photometric redshifts in 
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the range 2 < z < 3.5 have bhie optical to infrared col- 
ors, suggesting the presence of on going star formation, 
and that there were few, if any, candidates for red, non- 
star-forming, old galaxies at these redshifts. Could old 
objects with ~ 3 have faded past the NICMOS HDF-N 
detection limits? For a non-star-forming, passively evolv- 
ing galaxy at z — 2.7 with ages oi t = 0.5, 1.0, 2.0, and 
2.4 Gyr (the latter being the age of the Universe for the 
chosen cosmology), and assuming negligible dust extinc- 
tion, the stellar masses corresponding to the iJieo = 26.5 
limit of the NICMOS HDF N image arc M ~ 0.2, 0.5, 1.0, 
and 1.3 x 10^" -Mq, while the corresponding Vgoe — -H^ieo 
colors are 3.7, 6.1, 7.3, and 7.4 mag, respectively. Many of 
the LBGs in our sample have derived stellar masses that 
are greater than these limits. If old objects with compa- 
rable mass exist at these redshifts, then they should be 
detected in the deep, NICMOS data. Similarly, the lu- 
minosities of LBGs from our sample at 3 < ^; < 3.5 are 
large enough that even if their star formation were trun- 
cated at the observed epoch, they would still be detectable 
above the iJieo < 26.5 limiting magnitude ^ 1 Gyr later 
at 2.0 < 2 < 2.5. While these hypothetical old objects 
would be too red to meet the LEG selection criteria, we 
would expect photometric rcdshift techniques to identify 
them as high redshift, early-type-galaxy candidates. The 
"J-dropout" HDF-N-JDl might be such an object, but 



there are few other examples. In general, there is little 
difference in the color distribution between the spectro- 
scopically confirmed HDF-N galaxies with 2.0 < 2; < 3.5, 
and those with photometrically defined redshifts in the 
same range, and thus we conclude that there are few (if 
any) massive, old, non-star-forming objects in the HDF- 
N in this redshift interval. 

Although there are few candidates for red, non star- 
forming galaxies in this redshift range, there could still, in 
principle, be a significant mass component of old stars un- 
derlying the active star formation in the UV bright LBGs. 
The infrared photometry sets an upper limit to the allow- 
able total mass of the combined young and old popula- 
tions, and the fit to the full UV-to-optical SED sets a 
maximum mass ratio of old to young stars, or else the 
observed colors become too red. 

To set limits on the LEG mass contribution from 
old stars, we considered a simple, two-component star- 
formation history model. The "young" light from ongoing 
star formation is represented by the generic suite of models 
with monotonically decreasing or constant star formation 
histories from §4.1 (but restricted here to the Salpeter IMF 
and metallicities Z = 0.2 and 1.0 Zq). We then added a 
second component to the SED from an old stellar pop- 
ulation, formed in a (5-function burst that evolves pas- 
sively without further star formation. We assumed that 
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Fig. 19. — Schematic illustration showing the effect of adding a maximally old stellar population on the best— fitting SEDs of LBGs NIC 503 
(HDF 2-903.0) and 814 (HDF 3-93.0). For each object, the data points show the measured photometry with flux uncertainties and bandpass 
widths indicated by error bars. The solid black line shows the best-fitting "young" model SED for the Z = 0.2 Zq and Salpeter IMF models 
from Table 3. The red line shows the maximum allowable contribution to the SED from an additional, old stellar population whose age is 
that of the universe at the galaxy redshift. The mass ratios of the old-to-young stellar populations are given in the panels. The dashed black 
line shows the superposition of the two components. 



this component formed at z = oo, thus maximizing its age 
(1.75 — 3.15 Gyr for our adopted cosmofogy) and mass- 
to~hght ratio. The superposition of the two-component 
SEDs increases the galaxy luminosity predominantly for 
rest-frame wavelengths, Aq > 4000 A, as shown for two 
fiducial LBGs in Figure 19. We then refit each LBG with 
the full suite of two-component models using the same 
LBG photometry and uncertainties as in §4. This two 
component model lets us set an upper limit to the LBG 
stellar mass, using a stellar population with maximal M j L 
while still matching the observed UV-to-optical SED. 

In order to gauge the fit improvement due to adding a 
model component from an older stellar population, we ap- 
plied an i^-test on the minimum values of per degree 
of freedom derived for each galaxy using the two sets of 
models (the single-component and two-component star- 
formation histories). For most galaxies, the i^-test indi- 
cates that both models fit the data equally well. However, 
for a small fraction of the LBG sample, the addition of an 



old stellar component does improve the goodness-of-fit. 
Moreover, because the additional model constraint does 
not significantly degrade the fit to any of the LBGs in the 
sample, this may actually imply older stellar populations 
from previous star-formation episodes are indeed present 
(or at least for some objects). However, based on the lim- 
ited statistics, it is difficult to assign a significance to this 
hypothesis. One reasonable conclusion from this test is 
simply that there is not clear evidence for an absence of 
old stellar populations underlying the active star forma- 
tion in LBGs in this sample. 

The two-component models generally favor a slightly 
different range of parameters than those that best fit the 
single-age models. The best-fitting ages of the young stel- 
lar component from the two-component models are gen- 
erally somewhat smaller than those derived in §5.1, with 
correspondingly higher dust extinction. The mean val- 
ues for the LBG sample are t ^ 15 Myr (30 Myr), and 
Ai7oo « 2.4 mag (1.6 mag) for Z = 0.2 Zq (1.0 Zq). This 
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Fig. 20. — Constraints on the fraction of stellar mass in a hypothetical, old stellar component for two fiducial LBGs. The solid line and the 
hatched region show the best fitting value and 68% confidence region for the total stellar mass from the single component, "young" models 
with Salpcter IMF and Z = 0.2 Zq. The green shaded region shows the approximate 68% confidence region for the additional stellar mass 
allowable from a maximally old stellar population, in terms of the mass fraction = M{Old Stars)/ [A^ (Old) + (Young)] , as a function of the 
total stellar mass, A^(Old) + (Young), for best fitting, two component models. The vertical dashed line indicates the mass of a present-day 
L* galaxy. 



is simply a result of the "decomposition" of the observed 
SEDs into two components, where the "old" component 
influences the SED only past ~ 4400 A, with only a small 
amount of recently formed stars in the "young" compo- 
nent, which dominate the UV portion of the SED. 

The total stellar masses derived for these two- 
component models are larger than those from §5.1. This is 
due to the presence of the older stellar component, which 
has a large mass-to-light ratio and does not contribute 
substantially to the hght, but do contribute significantly 
more mass than do the younger stars. Figure 20 shows the 
relative mass fraction of old stars as a function of the total 
stellar mass (old + young stars) for the same two fiducial 
LBGs shown in Figure 19. For these galaxies, and indeed 
for most of the LBGs, there are viable models where the 
old stellar populations contain > 90% of the total stellar 
mass, while a minority population of young stars domi- 
nate the observed SED. However, for most galaxies, the 
total stellar masses (old and young stars) increase only by 
factors of a few over the estimates derived using the single- 
components models. In Figure 21, we compare the best- 
fitting stellar masses from one- and two-component star- 
formation history models. On average, the stellar masses 
from the two-component models are « 3x the masses de- 
rived from the single-component models. The (68%) up- 
per mass limits from the two-component models can be 
« 6 — 8x the best-fit single-component masses, with these 
old stars contributing < 20% of the light at rest-frame op- 
tical wavelengths. The fraction of mass that can be "hid- 
den" in an older stellar population increases with redshift 
as the observed bandpasses shift blueward, especially for 
z > 3 where the Balmer break shifts through the F160W 
bandpass. These "maximal A4/L," two-component star- 
formation history models set an upper bound to the to- 
tal stellar mass, and the results from the one- and two- 
component models bound the allowable range of LBG stel- 
lar masses for a given choice of IMF. Tighter constraints 



on the total stellar mass will require future observations 
at still longer rest -frame wavelengths. Deep observations 
with the SIRTF IRAC instrument will be able to measure 
the rest-frame if -band luminosities of LBGs, at these red- 
shifts, and should provide the best census of LBG stellar 
mass at these redshifts. 

5.3. Evolution of Lyman Break Galaxies 

The question remains, how do the LBGs fit into the his- 
tory of galaxy formation and evolution, and what is their 
fate at z = 0? As discussed in the introduction, several 
evolutionary scenarios have been proposed for LBGs. Re- 
gardless of the physical parameters that regulate star for- 
mation in these objects, these different scenarios predict 
substantially different end products for their evolution. 
The clustering observed for LBGs indicates that these 
galaxies are strongly biased relative to the overall mass 
distribution in the universe (Giavalisco et al. 1998, Adel- 
berger et al. 1998). This suggests that they are hosted by 
the most massive dark matter halos at these redshifts, and 
that they will ultimately become parts of fairly massive 
galaxies today. In this scenario, the LBGs may plausibly 
represent the formation stages of the bulges and spheroidal 
components of massive, present-day galaxies. In contrast, 
from their analyses of HDF-N LBGs, Sawicki & Yee (1998) 
concluded that the LBGs will form only a small fraction 
(~ 5%) of the stellar mass contained in a typical, L* galaxy 
today. That scenario is also consistent with one of the 
models suggested by Lowenthal et al. (1997), in which the 
LBGs could be small, sub-galactic clumps that will fade 
by the present day to represent present-day dwarf ellipti- 
cal/spheroidal galaxies. This question hinges on our un- 
derstanding of the fraction of stellar mass that formed has 
formed at these high redshifts. 

By fitting population synthesis models, we have found 
a range of star-formation histories that adequately match 
the observed photometry for each galaxy in our sample. 
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Fig. 21. — The best-fitting total stellar mass [= A4(01d Stars) + A^(Young Stars)] for each LBG, plotted against the mass from the 
single-component models. Each data point shows the range of total stellar mass from the two-component models as a function of the range 
of stellar masses from the single-component model fits. The error bars correspond to the 68% confidence range on each set of models. The 
solid line shows a one-to— one mass ratio, while the dashed line indicates the average mass ratio, X, of the best-fit values for the LBGs, where 
X = (2.9, 3.1) for the Z = (0.2, 1.0) Zq models, respectively. The dotted lines mark the stellar mass of a present-day L* galaxy. 



If we assume that models were to continue forming stars 
without change or interruption after the time at which 
the galaxies are observed, we may predict the subsequent 
evolution of these galaxies at lower redshifts. For each 
LBG, we selected all models within the 68% confidence 
region of the parameter space, i.e., those which best fit 
the photometric data. Because a range of different models 
provides acceptable fits to the photometry, there is there- 
fore a wide range of possible end-product masses. Star- 
formation models with t > t have already produced the 
bulk of their stellar populations by the redshift where the 
galaxy is observed. The colors and luminosities of these 
models at later times evolve essentially passively, and the 
stellar masses remain roughly constant, although some of 
the mass becomes locked up in stellar remnants or is re- 
turned to the interstellar medium as stars burn off the 
main sequence. Alternatively, models with t <^ t con- 
tinue to build stellar mass over their e-folding timescale, 
T, growing by factors of ^ r/t over their lifetimes. For 
most of the LBG sample, the best-fitting model remains 
significantly less massive than a present-day L* galaxy. 
However, it is also true that nearly every galaxy has possi- 
ble models (in the sense that they provide acceptable fits to 
the observed photometry) that will result in final stellar 
masses that are greater than those of L* galaxies today. 
Clearly, galaxies cannot continue to form stars continu- 
ously at high rates in this fashion for long, because it is 
unlikely that they have sufhciently large gas reservoirs to 
fuel such star formation indefinitely. 

Although the best-fit stellar masses for nearly all of 
the HDF-N LBGs are significantly smaller than that of 
a present-day L* galaxy, it is important to remember 
that most of these galaxies are also fainter than the "typ- 
ical" z w 3 LBGs studied in the ground-based survey 
of Steidel et al. This can be seen in Figures 17 and 18, 
where the masses are plotted against the UV luminosi- 
ties of the galaxies. The characteristic L* luminosity at 
1700 A for LBGs from Steidel et al. (1999) is approxi- 
mately Mi7oo = —21.0, which is brighter than most of the 



HDF-N LBGs. GiavaHsco & Dickinson (2001) have found 
that the fainter HDF-N LBGs cluster less strongly than 
do their brighter counterparts from the ground-based sur- 
veys, suggesting that they are generally associated with 
less massive dark matter halos. As a result, it might not 
be surprising if many of them are not in fact destined to 
become L > L* galaxies today. The stellar masses plotted 
in Figures 17 and 18 scale roughly linearly with the UV 
luminosity, and the intercept at M1700 = —21 is roughly 
10^°Alo , or ^ 1/lOth the mass of a present-day L* galaxy. 
As discussed in §5.2, we cannot discount the possibility 
that LBGs may have a component of older stars from pre- 
vious episodes of star formation. This could, in principle, 
significantly boost their stellar masses with little effect on 
their observed luminosities, but unless all LBGs have the 
maximum permissible amount of old stellar mass (i.e., un- 
less they actually formed most of their stars at extremely 
high redshift, z ^ 3), then it seems that LBGs with L* 
UV luminosities must generally have smaller stellar masses 
than their present-day L* descendents. 

For our adopted cosmology, the co-moving HDF- 
N volume between z = 2 and 3.5 is approximately 
27000/170"^ Mpc"^. From the if -band luminosity function 
of Cole et al. (2001), we would expect an average of 23 
galaxies with luminosities L > L* k within such a volume 
in the local universe. Modulo variations due to cluster- 
ing, their progenitors must be present within this redshift 
range in the HDF-N, at some earlier stage of evolution 
- perhaps in many smaller pieces, or even as yet uncol- 
lapsed. Of the galaxies that we have analyzed here, only 
a few have best-fitting star-formation histories that, if 
extrapolated forward without change to z = 0, would pro- 
duce all the stars found in an L* galaxy today. There 
are many other HDF-N LBGs without spectroscopic red- 
shifts (mostly fainter objects, see Figure 1) that we have 
not analyzed here, and some fraction of these may also 
grow to comparable masses at z = 0. In general, then, 
we can say that the best-fitting stellar population mod- 
els for HDF-N LBGs have masses smaller than those of 
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L* galaxies today, and that most LBGs are likely to re- 
main sub-L* without further episodes of star formation 
(perhaps fueled by gas infall) or substantial merging at 
lower redshifts. These HDF-N LBGs are less luminous, 
and thus arguably loss massive, than their counterparts 
from the larger ground based surveys, and it is unclear 
whether or not there are enough objects that could "pas- 
sively" evolve to produce the massive end of the present- 
day galaxy population without additional star formation 
or merging. However, these histories arc not definitive: for 
nearly every LEG, there exist allowable models capable of 
producing (or even overproducing) the total stellar mass 
of a present-day, L* galaxy. 

An interesting puzzle results from the fact that the ages 
of the 2.0 < z < 2.5 sub- sample of LEGs arc generally 
much shorter than the span of cosmic time across the red- 
shift range of our entire sample [t{2 < z < 3.5) ~ 1.5 Gyr]. 
Eecause the masses derived for LEGs at z ^ 3 — 3.5 are 
comparable to those at z ^ 2 — 2.5, we might expect to 
find a population of old, quiescent galaxies among the 
many HDF-N objects with (photometrically estimated) 
redshifts 2.0 < z < 2.5. Indeed, given the lifetimes we es- 
timate for the LEG star-formation episodes (with median 
stellar population ages t ~ 0.1 Gyr), such "old" galax- 
ies should outnumber the observed LBGs by a factor of 
^ 1.5 Gyr/0.1 Gyr = 15. We do not see such a population 
of old, high-redshift objects, given the narrow range colors 
observed for galaxies at these redshifts (sec Figure 1) to 
the flux limits of the NICMOS data. Moreover, Steidel et 
al. (1999) find that the UV luminosity functions of LEGs 
are indistinguishable (at least for luminosities around L*) 
at z ~ 3 and z ~ 4. The z ~ 4 objects must evolve into 
something, but again there is no evidence for a population 
of faded remnants at z < 3.5. Thus the fate of the the 
z ~ 3 LEG population remains somewhat enigmatic. 

One possible solution may be found if the simple star- 
formation prescription (i.e., constant or declining SFR 
with time) used in our models is incorrect for most objects. 
Galaxies might maintain a youthful appearance, despite 
long star-formation timescales, if the SFRs generally in- 
creased with time over this redshift range. This, however, 
may be difficult to reconcile with the apparent lack of evo- 
lution in the UV luminosity function between z = 3 and 4. 
Alternatively, multiple star-forming episodes with a duty 
cycle of < 1 Gyr might be the norm. In this way, the 
stars from high redshift LBGs are incorporated into their 
lower redshift counterparts, with intermittent episodes of 
star formation to rejuvenate the emergent UV-to-optical 
spectrum in order to avoid overproducing galaxies with red 
colors or old SED ages. Interestingly, when we fit LEGs 
using the two-component models (ongoing star formation 
plus a maximally-old stellar component, see §5.2), we fre- 
quently derive quite small ages for the "young" component 
(often t < 30 Myr), which may also be an indication that 
short bursts are the normal mode for star formation at 
these epochs. 

This complex, stochastic evolutionary history is remi- 
niscent of the collisional starburst models of Somerville et 

al. (2001), although in detail their model appears to pre- 
dict too broad a range of LEG stellar masses compared to 
the values we derive. In the general context of hierarchical 
models for galaxy formation, one might naturally expect 



that young galaxies are the sites of vigorous merging ac- 
tivity and gas accretion, which can trigger episodic (but 
frequent) star formation. The total stellar masses result- 
ing from these complex star-formation histories should be 
bounded at the minimum by the single component mod- 
els wc have used in this paper, and at the maximum by 
the two-component models with a maximally old stellar 
population. 

6. SUMMARY 

Using WFPC2 (C/3oo545oV'606/8i4), NICMOS 
(Jiio-ffieo), and ground-based iiTg-band data, we have 
investigated the stellar populations of 33 HDF-N LEGs 
with redshifts 2.0 < z < 3.5, reaching the following con- 
clusions: 

1. The LEGs have rest-frame UV-to-optical SEDs that 
are dissimilar to present-day, "Hubble Sequence" spiral 
and elliptical galaxies, but are comparable to those of local 
starburst galaxies with modest amounts of dust extinction. 
Examining complete galaxy samples based on photometric 
redshifts, we find little evidence for a substantial popula- 
tion of galaxies with redder colors at similar redshifts. 

2. By comparing the LEG data in seven broad-bands 
to stellar population synthesis models, we derive good 
constraints on the total stellar masses for the LEGs in 
our sample. The minimum value is derived using single- 
component star formation histories, with mean values 
that range from 3 (6) x 10^ Mq for Z = 0.2 Zq (1.0 Zq) 
and a Salpeter IMF to 7 x 10^ TW© using a Scalo IMF (for 
both low and solar metallicity models). For galaxies with 
L* UV luminosities, the inferred masses are ~ 10^'^ -Mq, 
or roughly 1/lOth that of a present day, L* galaxy. Fur- 
thermore, we constrain the maximum LEG masses by in- 
cluding a second model component, corresponding to a 
maximally-old (and hence maximum M / L) stellar popula- 
tion hidden in the glare of the ongoing star-formation. We 
find, on average, that a hypothetical old stellar population 
could in principle contain up to « 3 — 8 x the stellar mass 
of the young stars that dominate the observed SED. The 
stellar masses we derive are comparable to the simple kine- 
matic estimates for LEG masses from nebular-line widths, 
suggesting that those measurements underestimate the to- 
tal mass of the LEG dark matter halos. We find a fairly 
tight correlation between the derived stellar masses and 
UV luminosity, with a scatter of « 0.3 dex. 

3. The allowable (68% confidence) SED models gener- 
ally place only loose constraints on the stellar-population 
ages, dust attenuation, and prior star-formation history. 
For solar metallicity models, the best-fitting ages for in- 
dividual galaxies span a range of 30 Myr to ~ 1 Gyr, and 
have dust attenuation factors of — 3 mag at 1700 A. For 
models with sub-solar metallicity, however, a number of 
objects are well fit with very young ages (< 10 Myr) and 
high extinction (A1700 = 2.5 — 4 mag). The mean ages and 
extinctions for our LEG sample thus depend somewhat on 
the particular choice of IMF and metallicity, with ages of 
70 Myr (120 Myr) for 0.2 Zq (1.0 Zq) metallicity and 
Salpeter IMF, and 40 Myr (70 Myr) for 0.2 Zq (1.0 Zq) 
metallicity and Scalo IMF. The typical dust attenuation 
correction at 1700 A is a factor of 4.4 (3.0) for Z = 0.2 Zq 
(1.0 Zq), although the total, luminosity weighted mean 
flux correction for the sample as a whole is larger, 13.8 
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(4.6), due to larger fitted extinctions for some galaxies 
with the largest inferred SFRs. These values arc consistent 
with those from Steidel et al. (1999) based on rest-frame 
UV observations alone of nearly 1000 LBGs. These av- 
erage correction factors, while possibly applicable to the 
LEG population as a whole, are probably not very use- 
ful for individual galaxies due to the wide allowable range 
of parameter space for each LEG. We also find that no 
galaxy in our LEG sample is well fit by young, dust-frcc, 
"first-generation" stellar populations. Such objects must 
be quite rare at these redshifts, or too faint to meet the 
LEG color selection criteria, or must have SEDs that are 
significantly different from those produced by the stellar 
population models we use here. 

4. The extrapolated evolution of the best fit SED mod- 
els generally indicates that LBGs will remain sub-L*-sized 
galaxies today unless they undergo subsequent merging or 
gas accretion. However, for nearly all the LEGs, there are 
allowable (68% confidence) models capable of producing 
the mass of a present-day L* galaxy, and in some cases 
many times that of an L*-galaxy by ^ = 0. Based on 
the derived mass and age distributions, we suggest that 
star formation at high redshifts may be driven by dis- 
crete, recurrent events, which could result from mergers, 
tidal interactions, etc. Those starbursts must occur with 
a relatively short duty cycle (~ lO^Gyr/10^ Gyr), and a 
timescale between star-formation events of t < 1.5 Gyr, 
in order to keep rejuvenating the galaxy SEDs and ex- 



plain the apparent absence of non-star-forming, evolved, 
massive galaxies at z 2. The LEGs may indeed rep- 
resent the building sites of massive present-day galaxies, 
but the dominant stellar-mass assembly may occur via re- 
peating starburst events rather than quiescent, continuous 
star formation. 
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Table 1 



Observations of the HDF-N. 







Aoff 


Exposure Time 




PSF 


Camera 


Bandpass 


(/xm) 


(ksec) 


(AB)'' 


FWHM 


WFPC2 


F300W 


0.30 


153.7 


27.0 


0'.'14 




F450W 


0.46 


120.6 


27.9 






F606W 


0.60 


109.1 


28.2 






F814W 


0.80 


123.6 


27.6 




NIC3 


FllOW 


1.10 


12.6 


26.5 


0'.'22 




F160W 


1.60 


12.6 


26.5 




KPNO IRIM 




2.20 


82.4 


23.8 


1" 



^The WFPC2 and NICMOS limiting magnitudes are lOci in a 0.2 arcsec^ 
aperture. The IRIM JT^-band limiting magnitude is 50- in a 2 arcsec^ aper- 
ture. 
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Table 2 

HDF-N Lyman Break Galaxy Sample 



NIC ID"* 


WFPC2 ID'' 


z 


RA'^ 
(s) 


(mm:ss.s) 


HiQO 


Vqo6 — HiQQ 


Hl60 — I^s 


Ref. 


162^* 


2-585.2 


1.980 


49.74 


14:15.6 


24.46 


1.10 


< 1.02 


5,8 




2-585.1 


1.980 


49.82 


14:14.8 


22.44 


1.82 


0.32 


5,8 


no 


2-449.0 


2.005 


48.30 


14:16.6 


22.52 


1.21 


0.28 


1 


109 


2-454.0 


2.009 


48.23 


14:18.5 


23.49 


0.97 


0.08 


7 


1513 


3-875.0 


2.050 


37.07 


12.25.3 


23.06 


0.94 


-0.09 


8 


843 


3-118.0 


2.232 


54.71 


13:14.8 


23.56 


0.92 


0.19 


6 


503 


2-903.0 


2.233 


55.06 


13:47.1 


23.94 


0.70 


0.34 


2 


274 


2-525.0 


2.237 


50.09 


14:01.1 


23.58 


1.24 


0.51 


6 


67 


2-82.1 


2.267 


44.07 


14:10.1 


23.86 


0.67 


-0.12 


2 


561 


2-824.0 


2.419 


54.61 


13:41.4 


24.67 


0.63 


< -0.73 


2 


97* 


2-239.0 


2.427 


45.89 


14:12.1 


23.99 


0.75 


0.69 


6 


741 


2-591.2 


2.489 


53.18 


13:22.7 


23.99 


0.90 


0.29 


6 


989 


4-445.0 


2.500 


44.64 


12:27.4 


22.61 


1.55 


0.26 


3,8 


804 


4-639.0 


2.591 


41.72 


12:38.8 


24.53 


0.47 


0.52 


1 


1352 


4-497.0 


2.800 


45.74 


11:57.3 


25.11 


0.82 


0.32 


8 


782 


4-316.0 


2.801 


45.02 


12:51.1 


23.20 


1.17 


0.34 


8 


1357^^ 


4-555.1 


2.803 


45.41 


11:53.1 


23.21 


1.26 


0.28 


1 


1358^* 


4-555.2 


2.803 


45.30 


11:52.2 


23.49 


0.85 


0.45 


1 


522 


1-54.0 


2.929 


44.10 


13:10.8 


22.96 


1.44 


0.36 


6 


813^^ 


4-52.0 


2.931 


47.72 


12:55.8 


23.64 


1.29 


0.13 


2 


814^ 


3-93.0 


2.931 


47.86 


12:55.4 


23.65 


1.70 


0.40 


2 


J-UUO 




9 QfiQ 






9"^ 9fi 


n 1 s 


fl 7fi 


R 
u 


1541 


4-363.0 


2.980 


48.31 


11:45.8 


24.51 


0.88 


0.56 


2 


661 


2-643.0 


2.991 


53.42 


13:29.4 


24.15 


0.79 


0.02 


2 


273 


2-76.11 


3.160 


45.34 


13:47.0 


25.24 


0.24 


1.23 


2 


367 


2-565.0 


3.162 


51.19 


13:48.8 


24.15 


1.11 


-0.99 


6 


282 


2-901.0 


3.181 


53.58 


14:10.2 


23.97 


0.91 


-0.27 


2 


IIW'^ 


4-858.11 


3.220 


41.17 


12:02.9 


25.68 


0.53 


0.98 


1,2,4,8 


1115^* 


4-858.0 


3.220 


41.26 


12:03.0 


23.89 


0.77 


1.00 


1,2,4,8 


947 


3-243.0 


3.233 


49.80 


12:48.8 


25.51 


0.55 


1.64 


2 


284 


2-834.2 


3.367 


52.98 


14:08.5 


26.53 


0.65 


< 1.41 


8 


516 


2-637.0 


3.368 


52.74 


13:39.1 


24.88 


0.54 


0.38 


2 


553 


2-604.0 


3.430 


52.40 


13:37.8 


24.43 


0.80 


1.16 


2 



'^ID numbers refer to the catalog of Dickinson ct al. 2001. 
''ID numbers refer to the catalog of Williams ct al. 1996. 
<=RA add 12h36m; Declination add +62d (J2000) 

''Part of a close pair (separation < 1") split in our catalog and assumed (or confirmed) to have the 
same redshift. 

'^Multiple faint components in our catalog have been remerged (NICMOS IDs 96, 97, and 98). 

References. — (1) Steidel et al. 1996; (2) Lowcnthal ct al. 1997; (3) Cohen ct al. 1996; (4) Zcpf et 
al. 1997; (5) Elston et al., private communication; (6) Dickinson 1998; (7) Lowenthal, Simard, & Koo 
1998; (8) Cohen et al. 2000 
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Table 3 

Best Fit Stellar Population Synthesis Model Parameters 







Z 




Age 


T 


^1700^ 


log 


SFR 




NIC ID 


z 




IMF 


(Myr) 


(Myr) 


(mag) 


(Me) 


(TWq yr-i) 


Xo 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


162 


1.98 


0.2 


Salpeter 


9.1 


10 


3.1 


9.0 


68.9 


9.2 






1.0 


Salpeter 


90.5 


10 


0.5 


9.4 


0.035 


9.2 






U.z 


bcaio 


Q T 
O. ( 


lU 


O.i 


n A 


oni 
zUi 


y.i 






1.0 


bcaio 


80. 


20 


1 A 

1.0 


9.0 


4.12 


A 1 

9.1 


163 


1.98 


0.2 


Salpeter 


90.5 


10 


2.2 


10.5 


0.392 


9.5 






1.0 


Salpeter 


71.9 


10 


2.2 


10.4 


2.33 


9.7 






U.z 


bcalo 


1 no o 
iUz.U 


oU 


Z.4 


1 n 7 
iU. 1 


T^i n 
/ .U 


n 1 

y.i 






1.0 


Scalo 


d4.1 


1 A 

10 


o o 

2.2 


10.6 


7.56 


A A 

9.4 


110 


2.008 


0.2 


Salpeter 


37.0 


2000 


3.1 


10.1 


345 


3.3 






1.0 


Salpeter 


161.0 


2000 


2.4 


10.3 


148 


4.2 






U.z 


bcalo 


o4.U 


oU 


o n 

z.y 


1 n /) 
iU.4 


ooo 


O.O 






1.0 


bcaio 


OCA 

35.0 


1 A 

10 


2.2 


1 A /I 

1U.4 


AAA 

90.0 


A O 

4.3 


109 


2.008 


0.2 


CI 1 J. 

Salpeter 


7.6 


10 


3.1 


9.4 


217 


1.2 






1.0 


Salpeter 


55.0 


10 


1.2 


9.8 


2.75 


2.5 






U.z 


bcaio 


( .9 


oU 


O.i 


y.y 


Q1 n 


i.z 






1.0 


bcaio 


42.5 


1 A 

10 


1 A 

1.4 


1 A A 

lU.O 


15.1 


2.3 


1513 


2.05 


0.2 


baipctcr 


6.9 


10 


3.1 


9.6 


372 


1.8 






1.0 


Salpeter 


80.6 


10 


0.5 


10.0 


0.333 


5.8 






U.z 


bcalo 


D.y 


iU 


Q 1 

6.1 


inn 
iU.U 


iUo/ 


i. / 






1.0 


bcaio 


52.5 


1 A 

10 


1 A 

1.0 


1 A O 

1U.2 


7.71 


5.5 


843 


2.232 


0.2 


Salpctcr 


9.1 


10 


2.9 


9.4 


190 


1.1 






1.0 


Salpeter 


90.5 


30 


1.2 


9.9 


15.2 


0.6 






U.z 


bcalo 


oU.U 


Qn 

OU 


Z.4 


inn 
iU.U 


oi 1 
Zi ( 


n n 

u.y 






1.0 


bcalo 


50.0 


1 A 

10 


1 A 

1.0 


1 A A 

10.0 


7.55 


A C 

0.5 


503 


2.233 


0.2 


Salpctcr 


42.5 


10 


1.2 


9.5 


5.31 


3.2 






1.0 


Salpeter 


509.0 


00 


0.7 


9.8 


16.4 


4.6 






U.z 


bcalo 


tjy.U 


1 n 
lU 


i.z 


n Q 


io.4 


Q n 

o.y 






1.0 


bcalo 


203.0 


oo 


A T 

0.7 


A A 

9.9 


OA A 

39.9 


6.4 


274 


2.237 


0.2 


Salpeter 


47.5 


10 


2.2 


9.9 


7.60 


0.9 






1.0 


Salpeter 


641.0 


oo 


1.7 


10.2 


32.1 


2.1 






U.z 


bcalo 


4z.O 


1 n 
lU 


o o 
z.z 


1 n o 
iU.z 


on Q 
zU.O 


1 a. 
i.D 






1.0 


Scalo 


321.0 


00 


1.4 


10.3 


59.9 


4.0 


67 


2.267 


0.2 


Salpctcr 


143.0 


300 


1.2 


9.6 


25.3 


0.9 






1.0 


Salpeter 


404.0 


5000 


0.7 


9.8 


18.9 


1.0 






U.z 


bcaio 


1 no n 
iUz.U 


onnn 
zUUU 


1 o 
i.z 


n Q 


a. A a. 


n 7 

u. r 






1.0 


bcaio 


227.0 


3000 


0.5 


9.9 


32.4 


0.8 


561 


2.419 


0.2 


Salpctcr 


509.0 


OO 


0.5 


9.5 


7.31 


4.6 






1.0 


Salpctcr 


719.0 


1000 


0.0 


9.6 


4.59 


4.2 






U.z 


bcalo 


ZOD.U 


onnn 
oUUU 


n o 
U.z 


n 

y.o 


11/1 

ii.4 


A 7 

4. / 






1.0 


Scalo 


360.0 


00 


0.0 


9.6 


11.9 


5.2 


97 


2.427 


0.2 


Salpeter 


35.0 


10 


1.4 


9.5 


9.05 


6.8 






1.0 


Salpeter 


286.0 


oo 


1.0 


9.6 


16.9 


7.4 






0.2 


Scalo 


38.0 


20 


1.4 


9.7 


45.6 


7.2 






1.0 


Scalo 


128.0 


200 


0.7 


9.8 


34.0 


7.7 


741 


2.489 


0.2 


Salpeter 


57.1 


10 


1.4 


9.8 


2.16 


2.6 






1.0 


Salpctcr 


509.0 


7000 


1.4 


10.1 


26.3 


2.6 






0.2 


Scalo 


52.5 


10 


1.4 


10.0 


5.54 


2.7 






1.0 


Scalo 


286.0 


700 


1.0 


10.1 


34.4 


3.0 


989 


2.5 


0.2 


Salpeter 


2.5 


10 


5.8 


10.8 


23043 


7.9 






1.0 


Salpctcr 


4.8 


20 


5.1 


10.3 


3918 


8.1 






0.2 


Scalo 


2.6 


10 


5.3 


11.1 


42544 


9.6 






1.0 


Scalo 


4.2 


10 


5.1 


10.8 


12250 


8.6 
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Table 3 — Continued 







Z 




Age 




^1700^ 


log A4 


SFR 




NIC ID 


z 




IMF 


(Myr) 


(Myr) 


(mag) 


(Me) 


{Mq yr-i) 


Xo 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 




804 


2.591 


0.2 


Salpctcr 


3.6 


20 


2.9 


9.3 


505 


2.5 






1.0 


Salpctcr 


26.3 


10 


1.4 


9.3 


17.5 


1.4 






0.2 


Scalo 


3.5 


2000 


2.9 


9.8 


1795 


2.6 






1.0 


Scalo 


20.0 


10 


1.4 


9.6 


57.8 


1.5 


1352 


2.8 


0.2 


Salpeter 


102.0 


10 


0.0 


9.3 


0.009 


6.0 






1.0 


Salpctcr 


404.0 


200 


0.0 


9.5 


3.08 


5.2 






0.2 


Scalo 


161.0 


50 


0.0 


9.5 


2.77 


6.1 






1.0 


Scalo 


.■)()!). 


3()()() 


0.0 


9.7 


8.4() 


0.5 


782 


2.801 


0.2 


Salpctcr 


1.0 


10 


5.1 


10. G 


32518 


13.8 






1.0 


Salpctcr 


1.0 


10 


5.1 


10.4 


23486 


16.8 






0.2 


Scalo 


2.5 


10 


4.8 


10.9 


29417 


14.3 






1.0 


Scalo 


1.0 


10 


4.8 


10.8 


56860 


16.8 


1357 


2.803 


0.2 


Salpeter 


7.6 


10 


3.6 


9.9 


745 


0.4 






1.0 


Salpeter 


55.0 


10 


1.7 


10.3 


9.48 


0.6 






0.2 


Scalo 


7.9 


20 


3.6 


10.4 


2588 


0.5 






1.0 


Scalo 


57.1 


20 


1.9 


10.6 


114 


0.5 


1358 


2.803 


0.2 


Salpeter 


38.0 


10 


1.9 


10.1 


29.7 


1.5 






1.0 


Salpeter 


321.0 


7000 


1.4 


10.3 


72.1 


2.1 






0.2 


Scalo 


42.5 


20 


1.9 


10.3 


156 


2.1 






1.0 


Scalo 


143.0 


7000 


1.4 


10.3 


162 


2.9 


522 


2.929 


0.2 


Salpeter 


5.8 


200 


4.3 


10.3 


3095 


1.5 






1.0 


Salpctcr 


37.0 


10 


2.9 


10.5 


94.4 


2.1 






0.2 


Scalo 


5.5 


500 


4.3 


10.7 


9661 


1.5 






1.0 


Scalo 


38.0 


30 


3.1 


10.8 


806 


2.0 


813 


2.931 


0.2 


Salpeter 


102.0 


10 


1.0 


10.2 


0.062 


3.0 






1.0 


Salpeter 


143.0 


10 


0.0 


10.2 


0.001 


1.1 






0.2 


Scalo 


128.0 


10 


0.5 


10.3 


0.006 


2.5 






1.0 


Scalo 


161.0 


30 


0.0 


10.3 


3.11 


1.3 


814 


2.931 


0.2 


Salpeter 


10.0 


10 


4.3 


10.0 


569 


1.0 






1.0 


Salpeter 


80.6 


10 


1.9 


10.4 


0.882 


0.2 






0.2 


Scalo 


9.5 


10 


4.3 


10.4 


1565 


0.9 






1.0 


Scalo 


80.6 


20 


2.2 


10.6 


36.1 


0.2 


1063 


2.969 


0.2 


Salpeter 


50.0 


20 


0.7 


9.2 


7.26 


1.1 






1.0 


Salpeter 


453.0 


10000 


0.0 


9.4 


6.98 


1.1 






0.2 


Scalo 


40.0 


20 


0.7 


9.4 


20.5 


1.1 






1.0 


Scalo 


181.0 


oo 


0.0 


9.4 


15.5 


1.2 


1541 


2.98 


0.2 


Salpeter 


36.0 


20 


2.4 


9.7 


44.9 


0.3 






1.0 


Salpeter 


255.0 


10000 


1.7 


9.9 


33.2 


0.3 






0.2 


Scalo 


14.5 


10 


2.7 


9.9 


190 


0.4 






1.0 


Scalo 


114.0 


7000 


1.7 


10.0 


85.9 


0.4 


661 


2.991 


0.2 


Salpeter 


5.8 


10 


2.9 


9.5 


367 


1.2 






1.0 


Salpeter 


80.6 


70 


1.7 


9.8 


49.6 


1.6 






0.2 


Scalo 


6.0 


30 


2.9 


9.9 


1247 


1.2 






1.0 


Scalo 


71.9 


700 


1.7 


10.1 


163 


1.6 


273 


3.16 


0.2 


Salpctcr 


286.0 


10000 


0.5 


9.4 


8.67 


6.6 






1.0 


Salpeter 


641.0 


7000 


0.0 


9.6 


6.56 


6.1 






0.2 


Scalo 


143.0 


200 


0.2 


9.5 


14.7 


6.8 






1.0 


Scalo 


227.0 


oo 


0.0 


9.6 


14.2 


6.8 


367 


3.162 


0.2 


Salpctcr 


90.5 


10 


1.0 


10.0 


0.118 


10.5 






1.0 


Salpeter 


114.0 


20 


0.5 


9.9 


1.58 


10.4 






0.2 


Scalo 


128.0 


10 


0.2 


10.1 


0.004 


10.0 






1.0 


Scalo 


90.5 


10 


0.2 


10.1 


0.141 


10.1 
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Table 3 — Continued 







Z 




Age 


T 


^1700^ 


log 


SFR 




NIC ID 


z 




IMF 


(Myr) 


(Myr) 


(mag) 




(Mo yr-i) 




(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 






O 1 Q 1 

o.ioi 


U.z 


Salpctcr 


5.0 


on 
30 


O 1 


9. i 


/Oo 


4.7 






1.0 


Salpctcr 


128.0 


1000 


1.9 


10.1 


94.1 


5.2 






0.2 


Scalo 


5.5 


100 


3.1 


10.2 


2340 


4.7 






1.0 


Scalo 


64.1 


500 


1.9 


10.3 


255 


5.2 


1114 


3.226 


0.2 


Salpeter 


4.2 


10 


2.7 


9.0 


174 


3.9 






1.0 


Salpctcr 


161.0 


5000 


1.2 


9.3 


15.3 


4.1 






0.2 


Scalo 


4.0 


10 


2.7 


9.4 


542 


4.0 






1.0 


Scalo 


50.0 


00 


1.4 


9.4 


56.3 


4.2 


1115 


3.226 


0.2 


Salpeter 


453.0 


00 


1.4 
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^For starburst extinction law (Calzetti et al. 2000). 



